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The real-space renormalization group (RSRG) method introduced previously for the Brownian 
landscape is generalized to obtain the joint probability distribution of the subset of the important 
extrema at large scales of other one-dimensional landscapes. For a large class of models we give 
exact solutions obtained either by the use of constrained path-integrals in the continuum limit, or 
by solving the RSRG equations via an Ansatz which leads to the Liouville equation. We apply in 
particular our results to the toy model energy landscape, which consists in a quadratic potential 
plus a Brownian potential, which describes, among others, the energy of a single domain wall 
in a ID random field Ising model (RFIM) in the presence of a field gradient. The measure of the 
renormalized landscape is obtained explicitly in terms of Airy functions, and allows to study in details 
the Boltzmann equilibrium of a particle at low temperature as well as its non-equilibrium dynamics. 
For the equilibrium, we give results for the statistics of the absolute minimum which dominates 
at zero temperature, and for the configurations with nearly degenerate minima which govern the 
thermal fluctuations at very low-temperature. For the dynamics, we compute the distribution 
over samples of the equilibration time, or equivalently the distribution of the largest barrier in the 
system. We also study the properties of the rare configurations presenting an anomalously large 
equilibration time which govern the long-time dynamics. We compute the disorder averaged diffusion 
front, which interpolates between the Kesten distribution of the Sinai model at short rescaled time 
and the reaching of equilibrium at long rescaled time. Finally, the method allows to describe the full 
coarsening (i.e. many domain walls) of the ID RFIM in a field gradient as well as its equilibrium. 
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I. INTRODUCTION 



A. Overview 

Disordered systems |j],|2j are usually characterized by a complicated energy (or free energy) landscape, with many 
local minima representing metastable states, separated by saddles and barriers which govern the dynamics. In general, 
obtaining analytical information on the statistics of large scale extrema and saddles of a random energy landscape 
U{x) is an extremely hard task. In one dimension however, the situation is simpler since one has to deal with an 
alternating sequence of minima and maxima, and there is hope for analytical progress. Recently, the large scale 
extrema of the pure Brownian landscape and of the biased Brownian landscape have been studied with a powerful 
real space renormalization technique There the interest was in the non-equilibrium Arrhenius dynamics of a 

particle in such a landscape, the so-called Sinai model [||-^). At large time this dynamics is dominated by jumps 
between deep minima. The renormalization procedure thus consists in erasing iteratively the small barriers in the 
landscape, and retaining only the large scale, large energy features which control the large time dynamics. This 
procedure yields the asymptotically exact large time dynamics in the case of the Sinai model, and allows to obtain 
a host of exact results HQ. This real space renormalization technique has also been used in other one dimensional 
models, such as quantum spin chains Jll|Jl2]] , reaction diffusion with disorder ||lo|| and the random field Ising model 
, with the same fixed points associated to the pure or to the biased Brownian landscape. For the coarsening of 
pure ID $ 4 model at zero temperature from a random initial condition fL3| , another fixed point was found which 
corresponds to the successive elimination of smallest domains, the only variable there being the length. A question is 
then whether other energy landscapes can be studied using similar techniques. 



B. Definition and properties of the toy model 

A one-dimensional energy landscape of interest is the so-called "toy model" defined by the Hamiltonian 

Utoy(x) = |x 2 + V(x) (1) 

where the random potential V(x) is a Brownian motion with correlations 

(V(x) - V{y)f = 2\x - y\ (2) 

Here we have defined the model directly in the continuum limit, but the universality class contains all models with 
a small scale cutoff and the above behavior (^) at large scale. This model has been introduced by Villain et al. as 
a zero-dimensional toy-model for an interface in the random field Ising model [ fL4| . Let us briefly review previous 
studies of this model and explain some motivations to reexamine this model. 



1. Equilibrium : Imry-Ma argument and 'statistical tilt symmetry' 

For a single particle, the Boltzmann equilibrium problem at temperature T is defined by the partition function 
Z = j dxe- Vt °y^l T . Although it seems very simple at first sight, it is an interesting example of disordered system. 
Because of the quadratic well, the typical position of the particle is finite, but it fluctuates in a non trivial way from 
sample to sample. An Imry-Ma argument gives the following scaling for the typical scale of the absolute minimum 

x ~ (< x 2 > T =o) V2 = C/j- 2 / 3 (3) 

where we denote by the brackets < . > the thermal Boltzmann measure with potential (Q) and by an overbar the 
average over the random potential (0). The amplitude C will be computed here. Usual methods such as perturbation 
expansion (even summing all terms) ]T^] or iteration method ]l6[ ] completely fail to recover (^|) . A Larkin type pertur- 
bation [ pT{ with a small scale cutoff also predicts x ~ (1/m) instead of ([|). Within replica variational approximations 
of the present model (or within the solvable mean field limit of large embedding space dimension N — > oo, see below) 
this is cured at the price of spontaneous replica symmetry breaking (RSB) |2^,|3l|]. RSB is of course not expected 
to occur in the present exact model ([l] - |^) (see Ref. |Q for discussion of the cases where RSB may occur in a one 
dimensional model). 
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The toy model also possesses the 'statistical tilt symmetry' as other random held models jlq] or the directed polymer 
model [ ^9pC| ] , which imply exact remarquable identities between the disorder averages of the thermal cumulants of 
the position at any temperature, as a consequence from the statistical translation invariance of the random potential. 
These exact identities can be summarized by ]l8| 



In < e- Xx >= T— (4) 
2/x 



which is simply quadratic in A, so that the second cumulant reads 



T 



<x 2 >~<x> 2 = — (5) 
A* 

and the disorder-average of all thermal cumulants of order greater than two vanish! As emphasized in |L8| , in the low 
temperature regime, the result (^J) implies that the thermal fluctuations of the toy model are related to the presence 
of metastable states in rare disordered samples. 



2. N on- equilibrium dynamics : between Sinai and approach of equilibrium 



An Arrhcnius dynamics can be defined for the toy model, and provided /i is small so that barriers can be large, one 
can ask about universal long time properties. This was studied for the special case /i — where one recovers the Sinai 
model The dynamics for small but finite /j is interesting for several physical systems of the random held type, 

and has not been studied up to now, to our knowledge. Physically, one expects the diffusion to interpolate between 
a Sinai-type diffusion at short rescaled time and the approach of equilibrium at long rescaled time. 



C. Applications of the toy model 

1. Directed polymers in random media 

The toy model has also attracted a lot of attention in connection with the problem of D-dimensional manifolds 
in (N + D) random media, where it plays the role of a prototype model since it represents the simplest case D = 
and N = 1 pl] , p^ |. It has been moreover argued in |2^,|24| that the toy model actually describes the statistics of 
the end-point of a random polymer of length L in (1 + 1) with delta correlated noise in the limit L large with the 
correspondence [i = 1/L. It was found numerically in p5[ | that the correlation of the effective random potential for 
the directed polymer model behaves as (||) only at short distances but saturates at large distances. This analysis has 
been confirmed by a scaling argument in |l9| and was finally rigorously proven in [ p6| where the full statistics of the 
effective potential was exactly computed to be an "Airy process" , defined there. Thus the approximation of end-point 
of the directed polymer in (1 + 1) by the toy model is valid only locally. 



2. Burgers turbulence 

The toy model is related to the statistical properties of a velocity field v(x,t) evolving according to the Burgers 
equation 27 

dtv(x, t) + v(x, t)d x v(x, t) = i^d^v(x, t) (6) 
This equation can be solved via the Hopf-Cole transformation as 



1 / (x-vV 



f+°° d y e 2 "V 2t h0fe: 
U 0M) = : nx _ v) 2 ~ 3 (7) 

where 4>{y) = Jq dy'v(y', 0) is a Brownian motion in the case where the initial velocity field v(y' , t — 0) is a white 
noise process. In the inviscid limit v — > 0, the integrals are dominated by the saddle 
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which gives 



v(x,t) = X e(X;t) (9) 



The velocity distribution and the statistical properties of shocks (which appear at the points x where there are 
two degenerate minima in (|J)) are thus directly related to the properties of minima of the toy model (^) with the 
correspondence /i = l/t. 



3. Deterministic KPZ equation 
The toy model is also related to the relaxation of an interface evolving according to the deterministic KPZ equation 

us 

d t h(x,t) = vd*h+~{d x h{x,t)) 2 (10) 

from a random initial condition h(y,t = 0) which is a white noise in y. Indeed, the mapping ]29|]3C| 

v(x,t) = -\d x h(x,t) (11) 
transforms the deterministic KPZ equation (|l(i| ) into the burgers equation (^]), leading to 

h(x,t) = ^ ln / +00 ^ e -(^+-H (12) 



A 7-oo V&rvt 

Again in the limit v — » where the saddle method yields (^), the profile h(x, t) is made of parabolic segments joined 
together by cusps corresponding to the shocks of the Burgers equation (see for instance Fig.l in [p9[). 



4- Coarsening and equilibrium in some random field Ising model 

For definiteness, we will consider in the following an Ising spin chain, but one can study not only spin systems, but 
also for instance liquid helium meniscus on a rough substrate p3[ , which does also see random field disorder, and can 
in principle be probed in confined, quasi one dimensional geometry. 

Let us consider as in || the ID random field Ising system of Hamiltonian 

N-l N 

H = — J SfcSfc+i — ^ hkSk (13) 

fc=l k=l 

with spins Sk = ±1, and study its equilibrium and its low temperature Glauber dynamics. In the domain wall 
representation, two types of domains walls A, (+|— ), and B, (— 1+), in numbers Na and Nb, alternate at positions 
a a , 6 Q , respectively, and the Hamiltonian can be rewritten as 

N A N B 

H = H ref + 2J{N A + N B )-J2 U(a a ) + ^ U(b a ) (14) 

a— 1 a— 1 

with \Na ~ Nb\ < 1. The potential felt by domains walls of type B is 

k 

tf(*) = 2j>j (15) 

3=1 

while domain walls A feel potential —U(k). The dynamics proceeds by diffusion and annihilation of domains of 
opposite types. 
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In Q, the zero field hk = and the uniform applied field hk = H situations have been studied. The toy model 
corresponds to the other interesting case 

h^=kh (16) 



h 2 k - h k 2 = g (17) 

representing a field gradient with zero net applied field (such gradients are easy to produce experimentally) . Indeed, 
it is easy to see that defining the scaled length x — 2gk, each domain wall B sees the toy model energy landscape 
U(k) — U toy (x) with U toy {x) defined in (jlj) in terms of a unit Brownian and 

*-£ (18) 

For a positive field gradient (h > 0), the domain walls of type B (— 1+) will be attracted towards the minima of 
U t0 y(x), and will thus be confined in the neighborhood of the point x = where the applied field vanishes, whereas 
the domain walls of type A (+|— ) will be attracted by the local maxima of Ut oy (x) or will be expelled towards the 
boundaries. 

The equilibrium of the toy model considered above and its dynamics involving only one particle corresponds to 
a single domain wall B. If the field gradient is ft. > and if the system is initially prepared with the "natural" 
boundary conditions, namely spin far on the left are (— ), far on the right (+) and with a single domain wall of type 
B in between, it remains such, until equilibration time Tlnt eq o = 2 J (or Thit eq Q — J if domains can be created at 
boundaries). During that window of time the domain wall diffuses in the toy model landscape and reaches thermal 
equilibrium at time t eq , a time scale which can be much smaller than t eq o if J is sufficiently large. Another interesting 
situation corresponds to unstable single domain wall dynamics and is obtained if the field gradient is applied against 
these boundary conditions, e.g. a single A domain initially. The domain is then expelled towards either side after a 
time t exp which again can be well before t eq Q if J is sufficiently large. 

Finally, the full, many domain coarsening can be studied and is expected to exhibit a large degree of universality 
in the small h, large Imry Ma length 4J 2 /g, small temperature T <C J, large time limit as in [Q. It is interesting 
to study because of the competition between the disorder which tends to slow the coarsening and the field gradient 
which tend to accelerate the process. 



D. Goal and Results 

To summarize, the toy model (Q) is thus interesting from several points of view. For the statics, it is directly related 
to more complex problems ( directed polymer, burgers turbulence and KPZ equation) in some particular limits, as 
described above; on the other hand, it represents in itself a 'toy' disordered system that presents many properties of 
more complex models, such as the presence of metastable states and the identities of the statistical tilt symmetry. 

Up to now, the landscape U (x) of the toy model has been characterized by two kinds of approach. There has been 
on one side some studies on the absolute minimum of the toy model, either in probability theory |3^] and in the 
context of the Burgers equation j2£| , and on the other side a study |52[ on the statistics of all extrema of the model 
, but that does not distinguish between the important ones (the deep ones) and the large number of unimportant 
small scale extrema, nor does it yield any information on joint distributions. None of these studies addresses the 
question of the joint distributions in energy and position of the set of all deep extrema, neither the one relevant for 
the statics (several minima) and even less the one relevant for the dynamics (maxima, minima and barriers). 

The goal of this paper is to show that the renormalization method can be generalized to solve a broader class of 
one dimensional landscape than the previously solved Brownian landscape jl2|,|} and that it is of deeper fundamental 
interest as an alternative to more conventional probability theory techniques. The class explicitly solved here includes 
the toy model ([!]) landscape. We formulate the RSRG for a general one-dimensional landscape as an exact method to 
obtain the joint probability distribution, in total number, position and value of energy on the subset of the important 
deep extrema. More precisely, the renormalized landscape at scale T is defined by the subset of extrema of the initial 
landscape which have survived after the erasing of all barriers smaller than T. This is equivalent to require that in 
a segment between a maximum U (M) and the next minimum U (m) of the renormalized landscape, any two points 
(x,y) with M < x < y < m must satisfy the constraint that the energy difference U(y) — U(x) is smaller than T 
(and the symmetrical condition for a segment between a minimum and the next maximum). The joint probability 
distribution of this subset of important extrema contains a tremendous amount of information, including in particular 
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the distribution of absolute extrema, of secondary extrema, of largest barriers etc... Of course the general case is 
intractable analytically, but our aim is to identify solvable models and universality classes, given by some continuum 
limits. 

To be more specific, we will consider the continuum models where the potential U(x) is a stochastic process defined 
by the Langevin equation 

^M=F[U(x),x}+ V (x) (19) 
ax 

rj(x) being a delta correlated white noise. The toy model (|l]) corresponds to the case where the force doesn't depend 
on U (x) and is linear in x 

F[U(x),x] = F [x] = fix (20) 

Another cases of interest are the stationary landscapes, where the force in the Langevin equation ( jigj ) is independent 
of the space x 

The case F[U] — is the pure Brownian landscape, and the case F[U] = F > is the biased Brownian landscape. 
These are the only two cases where the RSRG has been applied up to now [fi2p|], but we will show here that it 
generalizes to include any F[U}. Note that solving the landscape problem here also amounts to solve the problem 
of the joint distribution of extremal points in the trajectory of a particle in a one dimensional force field F where 
x is interpreted as "time" and U as the position of the particle. For all cases the Markovian character of the 
Langevin equation implies that the landscape measure can be written as product of blocks measures which satisfy 
closed real space renormalization equations. For the specific cases of the toy model ( p0| ) and of stationary landscapes 
(pl|), we obtain here explicit expressions for the blocks measures. This is achieved by two complementary methods 
very different in spirit. On one hand, for the above continuum models we are able to directly compute the blocks 
using some constrained path integrals. On the other hand, we can also start from the RSRG equations for the blocks 
and assume some Ansatz which solves these equations and obtain, interestingly via the Liouville equation, analytical 
forms for the blocks. As is customary in searching for RG fixed points, this can be done without specifying the initial 
landscape, and as a result the solutions look a priori more general. We show that these do contain the universality 
classes corresponding to (|l9|), with a result which coincide exactly with one obtained by the direct path integral 
method. Our study thus makes more precise the interesting relations between the powerful RSRG method and more 
conventional probabilistic approach using stochastic calculus. 

For the case of the toy model, the full solution for the joint probability measure of the renormalized landscape 
is obtained in this way explicitly in terms of Airy functions. We then use this exact measure of the renormalized 
landscape to study in details the equilibrium and the non-equilibrium dynamics of the toy model (|l|) and give some 
results for the RFIM in a field gradient. 

For the equilibrium, we study the statistical properties of the absolute minimum which dominates at zero temper- 
ature. We obtain simple exact expressions for the moments of the position of a single particle at equilibrium and and 
compare with the results of approximate methods, namely the gaussian variational replica study p^j3l[ | and of the 
vector breaking of replica symmetry |38[ . We also compute the properties of configurations with nearly degenerate 
minima which govern the thermal fluctuations at very low-temperature, and we prove that the rare configurations pre- 
senting two nearly degenerate minima do not only give a finite contribution to (||) , but actually are entirely responsible 
of the full result (|). 

For the dynamics of a single particle in the toy energy landscape (e.g. a single RFIM domain wall), we compute 
the distribution over samples of the equilibration time t eq or equivalently the distribution of the largest barrier 
T eq = T In t eq in the system. Our results are expressed in terms of the dimensionless scaling parameter 



leq 



^) 1/3 Tlnt eq (22) 



and we give in particular the two asymptotic behaviors of the distribution for ~/ eq small and large. We characterize 
the rare configurations presenting an anomalously large equilibration time which govern the long-time dynamics. We 
also compute local properties of the renormalized landscape, such as the probability density of a minimum at scale 
r = Tint, and the the disorder averaged diffusion front for a single particle, which interpolates between the Kesten 
distribution of the Sinai model at short rescaled time and the reaching of equilibrium at long rescaled time. 



G 



These quantities can be translated into results for observables of the RFIM in a field gradient, assuming convergence 
to the full state where each maximum of the renormalized landscape contains an A domain and each minimum a B 
domain. This convergence was proved in || and also holds here in range of parameters of interest (small (j,, large 
time t) and with a random initial condition. As time r = Tint increases, decimations occur corresponding to pairs of 
A — B domains annihilating. Thus initially, many domain coarsening take place. Depending on the scaling parameter: 

1/3 



leqO = 2J (- 



(23) 



the equilibrium state, reached at time Thit eq o = 2 J will contain either a single domain wall B (for 7 eg o — ► +oo) or 
several additional pairs A — B (infinitely many as 7 eg o — > 0). The behaviors of the sample to sample probabilities 
that a single domain wall B remains is obtained. A formal solution for the total magnetization is easily written, and 
an explicit calculation is performed for large 7 eq o where a single domain wall remains. Indeed in that limit the total 
magnetization is simply M — 2x R where x R is the absolute minimum of the toy landscape. 

The paper is organized as follows. The Section [n] is devoted to the general definitions and properties of the 
renormalization procedure for an arbitrary one-dimensional landscape. In Section [II, we compute the probability 
measure of the renormalized landscape for a large class of models. We then study in details the toy model : we 
construct the probability measure for the renormalized landscape in Section IV, we study its equilibrium properties 
in Section 0, and the non-equilibrium dy nam ics in Section VI. Some results for the RFIM in a field gradient are 
also discussed in Section Jv|. The Section VII contains our results for the case of stationary landscapes. Finally, to 
facilitate reading of the paper, we give the more detailed definitions and derivations in the several appendices. 



II. RENORMALIZATION PROCEDURE FOR AN ARBITRARY ID LANDSCAPE 

This section contains the precise definitions of the measure for the initial landscape and for the renormalized 
landscape. The probability measure of the renormalized landscape satisfies by construction a set of RG equations for 
the joint distribution of extrema. 



A. Initial landscape 

The general model studied in this paper is the following. We consider one dimensional random landscapes U(x) 
which for definiteness in this section are continuous in the energy U , discrete in space x and defined on a finite 
interval x = 1, 2, ...N. A realization of the random landscape is thus simply given by a set of random variables U x , 

x = 1, 2, ...TV and the original problem is defined by some probability measure noted 

V[U]DU = V{Ux,...U N ) dUi...dU N (24) 
The normalization with the fixed boundary conditions U(x\ = Xl) = Ul and U{xm = xr) = U R reads 

Z(U L ,x L ;U R ,x R ) = Jv(U L ,U 2 ,...U N . 1 ,U R ) dU 2 ...dU N _i (25) 

We are interested in the statistics of the extrema of this landscape. For definiteness, we choose here the convention 
that there are infinite potential outside the interval, i.e. U(x = 1~) = U(x = N + ) = +oo, and that N — 2K — 1 
is odd. We then define the probability iVfij/ (Ul] h>i, X\; . . .u-^n-ii x 2n —x\ x R , Ur) that the initial landscape has n 
local minima situated at positions Xx, Xs-X 2n -i with the potentials u\, ...u 2n -i and (n — 1) local maxima situated at 
positions X2, X4 : ..X2 n ~2 with potentials 112, ...ii2n-2- They are of course alternated x\ < X2-. < X2 n -i and the values 

xi = xl and X2 n -i = X R = 2-ftT — 1 are allowed. An explicit and more precise definition of these objects is given in 
Appendix with a derivation of their normalization property 

K 

Z(U l ,x l ;Ur,x r ) =2J / du 1 du2...du 2n -i6(U L - u 1 )6{U R - u 2n -i) 0(u 2 i - u 2 i-i)0(u 2 i - u 2t+ i) 

n—l i— l,n— 1 

^ (Ul,x l ;ui,xi; . . .u 2n -i,x 2n -i; U R ,x R ) (26) 

x L <x 1 <x 2 <x 3 <x 2n -i<Xu 

The objects are usually extremely complicated and contain a lot of irrelevant information. To single out the 

important, deep extrema we now define the renormalized landscape where all small barriers have been eliminated. 
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B. Construction of the renormalized landscape 



The renormalized landscape at scale T is defined by the list of positions and energies of a subset Sr of extrema of 
the initial landscape. For r = 0, So is the set of all extrema of the initial landscape. For T > it is the subset of 
extrema that survive when all barriers smaller than T in the system have been erased. This is shown in Fig 1. This 
subset Sr is again made of alternating minima and maxima. It is convenient to call descending bonds the segments 
between a maximum and the next minimum both in Sr (and respectively ascending bonds the segments between a 
minimum and the next maximum). The constraint obeyed by the subset is that any two points on a bond [x,y] 
with x < i < j < y must satisfy that their energy difference Uj — Ui is smaller than T if the bond is descending, and 
symmetrically, the energy difference Ui — Uj must be smaller than T if the bond is ascending. 



U(x) 




f= 1 5 10 15 20 X=25 



U(x) 




\ X l X 2 X 3 X 4 X 5 X R 

FIG. 1. top One realization of the initial random landscape U(x) drawn from the probability distribution V[U] with N = 25. 

It contributes to N-ffl (xl, ul; ui,xi; . . . U2n-i,:E2n-i; xr, ur) with n — 6 and (asi, ..in) = (2, 3, 4, 6, 10, 11, 12, 17, 21, 22, 24) 
since only 2n — 1 = 11 points of the initial landscape are actually extrema. bottom The corresponding realization 
of the renormalized landscape where all barriers smaller than some given V have been removed. It contributes to 
iVp 2n ' (xl,ul',Ui,xi; . . . u 2n >_i, x 2n >_i; xr, ur) with n! = 3. As in the top figure, the straight lines are drawn only for conve- 
nience, the landscape being defined solely by the set of (xi, m), i = 1, .., 2n' — 1 = 5 of remaining extrema, together with the 
right and left boundary points. 

The main object of interest is the probability measure of the renormalized landscape. We consider the probability 
Nr (xl, ul; u±, x±; . . . U2n-i, %2n-i] xr, ur) that the renormalized landscape has n local minima situated at positions 
Xt, x$..X2n-i with the potentials ui,...U2n—i an d {n— 1) local maxima situated at positions xi,x^..xm~i with po- 
tentials U2, ...«2 n -2- Again, more precise and explicit definitions are given in Appendix |X[ with a derivation of the 
normalization 

K 

Z(Ul,xl;Ur,xr) = 22 / dui...du 2n -iO(UL - u\)0{Ur - u 2n -i) 0(u2i - U2i-i - T)6(u2i - «2»+i - T) 

n— 1 i— l,n— 1 

^ N r 2n) (U L ,x L ; ui,xi; . . . u 2n -i, x 2n -i; Ur,x r ) (27) 

I[<I1<I 2 <I3<I2„-1<IB 



The normalization of the measure for the renormalized landscape thus coincides with the normalization of the measure 
for the initial landscape ( p5| , p6| ). 

The N r 2n ^ satisfy the RG equation (see Appendix |A]) 



drN T 2n) {x L ,u L ;ui,xi; . . .u 2n -i,X2n-i;x R ,UR) 
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xr,ur) 



/ dy 1 dy 2 du' Nf n+2) (x L , u L ; u' , y x ; u' + T, y 2 ; u 1; x x ; . . . u 2n -i, x 2n -\\ %R, ur) 

>x L Jyi J Ui 

rxai+i rXM+i rU2i—T 

f ^ / dy x I dy 2 du'Nf n+2) (x L , u L ; . . . ; u 2ll x 2l ; u', y x ; u' + T, y 2 ; u 2i+1 , x 2l+1 ; 

J X2i J yi Jtl2i + 1 

f / ^2/1 / ^2/2/ du'N^ n+2) (x L ,u L ;...;u 2 i-i,x 2 i-i;u' + T, yi;u' ,y 2l u 2i ,x 2i ; . . .^xr.ur) 

i=l J x 2i-l Jyi JU2i-l 

x R rx R r+oo 

dyi dy 2 du 1 'Nf n+2) (x L ,u L ;ux,Xx; . . . u 2 „_i, a; 2n _i; u' + T,yi;u' \y 2 ;x R ,u R ) (28) 



where, for convenience and future applications to continuum limits, we have used the integral symbol instead of the 
sum, i.e. dz instead of X)z=y+i- These RSRG equation describe the decimation upon increasing V and are thus 
a generalization to an arbitrary landscape of the standard RSRG equations considered in the case of a Brownian 
landscape (see the formulae (197) and the figure 8 in Ref 0). It simply expresses that as T is increased, pairs of 
extrema separated by an energy barrier V are eliminated. Each term in ( p8|) corresponds to possible positions of this 
pair, as represented in Fig 0, the first and last one corresponding to edge bonds. 

The aim of this paper is to find solvable cases where analytical solutions of (28) are possible. The general case with 
arbitrary correlations in the original landscape can only be studied by numerical implementations of the decimation 
procedure. This has been done, for instance, in the case of a logarithmically correlated landscape in j35j. However, 
for the initial landscapes that are defined by an initial local measure, we expect that the measure for the renormalized 
landscape can be decomposed as a product of blocks associated to bonds as we now describe. 




FIG. 2. Terms in the RSRG equation ( Pq ) corresponding to the elimination of two local extrema separated by a potential 
difference V, for a descending bond (top) and ascending (bottom) respectively 



C. Markovian landscapes and Block-product measures 



From now on, we will consider "Markovian landscapes" , corresponding to the cases where the probability measure 
of the initial landscape (M) is a product of block measures of the form: 



N-l 



V[U] = MUi)*n(U n ) II ex.p(-Si[U i+ uUi\) 



(29) 



It follows immediately that the probabilities of the renormalized landscape can be written as a product of block 
probabilities which depend only on the endpoints of the bonds block, all the information inside the bonds have been 
in a sense integrated out ( A15) : 
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N^ 2n) {xl,uc;ui,xi; . . . u 2n -i,x 2n -i;x R ,u R ) = E£(x L ,u L ;ui,Xi)B T (ui,x 1 ;u 2 ,x 2 ) 
B+(u 2 ,x 2 ;u 3 ,x 3 ) . . .B+(u 2n - 2 ,x 2n - 2 ;u 2n -i,x 2n -i)E~ {u 2n -i,x 2n -i;x R ,u R ) (30) 

Note that the blocks B v (u 2l x 2 ; 113, X3) themselves do not in general satisfy any normalization. The only constraint 
is the global normalization: 

K 

Z{u L ,x L ;u R ,x R ) = Y^TriE+iB-B+y^E-) (31) 

71=1 

where we have used an obvious matrix notation with respect to (x, u). Each term in this sum, 
Tr(E + (B~ B + ) n ^ 1 E^), represents the probability that the renormalized landscape has exactly (2n — 1) extrema. 

The insertion of the decoupled form ([5(j|) into the RSRG equation ( p8| ) yields the following RG equations for the 
block measures 

d r E^(x L ,u L ;u,x) 

rx z'+oo 

dyi dy 2 du'E^(xL,U L ;u' ,yi)B^(u' ,yi;u' + T,y 2 )B^(u' +T,y 2 ;u,x) (32) 



II ~.H\ 



drB T [u, x; u , x' 

c-x pu — r 

dyi dy 2 du'B^ (u,x;u' + T,y 1 )B^ (u' + T,y 1 ;u' ,y 2 )B^ (u 1 ,y 2 ;u" ,x") (33) 

J Vi 

drB+(u,x;u",x") 

fX rU — T 

dyi dy 2 du'B+(u,x;u', yi )B^(u', yi ;u' +T,y 2 )B+(u' +T,y 2 ;u",x") (34) 



in 



drE r (u, x; x R , u R ) 

C~\~OQ 



rx fx r-j-00 

/ dyi dy 2 dv! ' B^ (u,x;vf +T,yi)B^ (u' + T,y 1 ;u' ,y 2 )E^ (u' ,y 2 ;x R ,u R ) 
J xl J yi ** u 



(35) 



Again note that the factorized form and these RG equations hold either for discrete space (with the abovementioned 
correspondence between the integral symbol and the sum) as well as in the continuum limit where the space cutoff 
goes to zero. 



III. EXPLICIT SOLUTIONS FOR A CLASS OF MODELS 
A. Soluble models 

To be more specific we now consider the continuum models of Markovian landscapes where the potential U(x) is 
a stochastic process defined by the Langevin equation (|l9|). This process corresponds to functional integration with 
the Fokker-Planck measure 

T[U}~eM- I dx(\(^l - F[U(x),x}) 2 + ld u F[U{x),x}) (36) 
J 4 dx 2 

In the following, we will give an explicit solution for the landscape renormalization for all the cases where the 
Fokker-Planck measure can be rewritten as a Schrodinger measure 

V[U] ~ exp(- J dx{\{^± f + V[U{x)}) (37) 

up to boundaries terms. The two important cases are the toy model defined in (pi) and the stationary landscapes (|2 
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1. Toy model 

The toy model (Q) corresponds to the case where the force in the Langevin equation (19) is given by (p0|). The 
Fokker-Planck measure for the toy model can be be rewritten as a Schrodinger measure up to boundary terms as 
follows 



f X2 1 dlJ P X2 1 dli 1 

cxp(- / dx(-(^- - fix) 2 ) = e h ("i.*i)- h ("*.*2)exp(- / dx{-(-^—) 2 H — fJ-U)) (38) 
' J X1 4 dx 2 



4 dx 



where the boundary terms reads 



2 

h(u, x) — — — xu + — x 3 (39) 



The associated Schrodinger potential is linear in U 



V[U] = ~fJ,U (40) 



2. Stationary landscapes 

The stationary landscapes correspond to the cases where the force in the Langevin equation ([l9]) is independent 
of space ( |2l| ) . The transformation from Fokker-Planck measure to Schrodinger measure up to boundary terms is as 
follows: 



f X2 1 dJJ 1 
cxp(- £ dx(-{— - F[U]f + -F'[U])) 



e Uw[u( Xl )]-w[u( X2 )]) gxp( _ f X2 dx (-& 2 + -F[U] 2 + -F'[U])) (41) 

Jxi ^ ^ ^ 



The associated Schrodinger potential reads 



V[U] = \f[U] 2 + l -F'[U] = \W[U] 2 \w"[U] (42) 



B. Properties of the blocks for a Schrodinger measure 

For a Schrodinger measure ( |37|) , the important simplification comes from the translation invariance of the blocks 

B T (it, x; u" , x") = _Bp (u; u"; x" — x) (43) 
E£(u Ll x l ; u,x) = E£ (ul, u,x - x L ) (44) 
E~(u,x;u R ,x R ) = E+(u,u R ,x R - x) (45) 

which will allow us to obtain explicit expressions for the Laplace transforms 

B^{u;u";p)= / dxB^ (u; u"; x) (46) 

Ey{u;u";p)= / dxE r (u; u"; x) 
Jo 

The symmetry x — > — x moreover implies: 

B+(u,x;u",x") = B-{u,-x;u",-x") (47) 

so that in Laplace we have the simple relations 

B+{u-u"-p)=B-{u"-u-p) (48) 
E+(u;u";p) = Er(u";u;p) (49) 

To determine these functions, we now present two methods which are rather different in spirit and complementary. 
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C. Explicit solution via a path-integral approach 



For the Schrodinger measure, the block B r can be represented by a path-integral 



B^(u ,x ;u,x) = DUex P (- d y {-( — Y + V[U])Q^{x 0l x;U) (50) 

JU X0 =u„ Jxo 4 d V 

which is defined for u > uq + T and where the symbol 8p(so,s; ?7) constraints the paths to remain in the interval 



]uq, u\ for xq < y < x and to not perform any returns of more than T. As shown in Appendix (C 1), this constrained 
path-integral can be expressed in Laplace transform with respect to (x — xq) as 

B~(u ,u,p) = dxe~ px B^(u ,x ;u,x)= — — — -exp(-/ du'd 2 InK (v! - T, u',p)) (51) 

Jx K{uo,uo + T lP ) J un+r 



K(u ,u + T,p) J Uo+r 
i ru-T 



K(u — r, u,p) 



exp(/ dvd x \nK{v,v + T,p)) (52) 



no 



(the two expressions are equivalent via an integration by part) in terms of the function 

K(u,v,p) = -^r^((j)x{u,p)4>2{v,p) - <j)2{u,p)4> x {v,p)) (53) 

where (f>i(u,p) and <f>2(u,p) are two linearly independent solutions of the equation 

d 2 u F - V{u)F = pF (54) 

of wronskian ( independent of it ) 

w{p) = (j)i(u,p)(j)' 2 {u,p) - 2 (u,p)<?!>i(u,p) (55) 



Similarly for the edges, the final result reads ( see Appendix C 1 



E T (u,ur,p) = exp( / d% hnK(u',u' + T,p)du') (56) 



D. Explicit solution from the RSRG equations 

The RG equations for the blocks B (|35| ) are simpler in Laplace variables : 

d v B- (u;u";p) = \ du'B^(u;u' + T;p)B+(u' + T;u';p)B-(u';u";p) (57) 



and the edges satisfy 

r+oo 



(58) 



nil— I 

drB+(u;u";p) = / du'B^ (u; u'; p)B^ (vf; u' + T; p)B$ (vf + T; u"; p) 

J u" 

Using the symmetry (|9|), we obtain that B^ (u"; u;p) satisfies the closed RG equation 

nu"-V 

d r B^{u;u";p) = du'B^(u; v! + T;p)B ] : (u 1 ; u' + T;p)B^(u'; u";p) (59) 

J u 



r-t-oo 

dpE^ (u,ur,p) = / du' B^ (u;u' + T;p)Bf (u 1 ;u' +T;p)E^ (u' ,Ufi,p) (60) 

J u 



As shown in Appendix (C 2), a class of solutions of the RSRG equation ( |59| ) can be found with the factorized Ansatz 

B- (u;u";p) = A$(u,p)A*(u" ,p) (61) 
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provided that the function $ defined by 



d 1 d 2 ^{u,u + Y)=B^{u-u + T;pf (62) 
satisfies the following Liouville equation 

<9i<9 2 $(ui, u 2 ) = exp(-2$(u!, u 2 )) (63) 



The general solution, discussed in more details in Appendix (C2), can be written 



$(ui, u 2 ) = In 



(64) 



(1 2) (1 2) 

where 4> L ' (u) and ip R ' (u) are respectively two pairs of two linearly independent solutions of the two Schrodinger 
equations: 

*Pl( u ) = V L (u,p)ip L (u) (65) 
^{u) = V R {u,p)^ R {u) (66) 

where Vl(u,p) and V R (u,p) are arbitrary potentials, and where the constants wl,r are the two wronskian, i.e wl.r = 

The results of the path-integral approach (^Cj) are recovered with the following choice: 

V L (u,p)=V R (u,p)=V(u)+p (67) 
with the identifications ip^ = ip^ , — ip^ 1 and 

$(iH,u 2 ) =\tlK{u 1 ,u 2 ) (68) 

Thus by solving the RSRG equations with the following properties: (i) factorization in blocks (ii) statistical trans- 
lational invariance and parity of the blocks (iii) an additional factorization property in Laplace variable (|6l]) , we have 
found general solutions involving two arbitrary potentials Vl(u,p) and V R (u,p). This method does not use the knowl- 
edge of the initial landscape they correspond to. This is a familiar feature of RG equations, fixed points can be found 
under symmetry assumptions but finding the basin of attraction is harder. We have showed that the special choice ( |6^ ) 
corresponds to the case of initial Schrodinger measure ( |50| ) that we have independently solved via the path-integral 
method. It would be of course of great interest to understand whether other choices for the potentials Vl(u,p) and 
V R (u,p) yields sensible solutions for the landscape problem. However this analysis goes beyond the scope of this paper 
and it is left as an open question for future investigation. Note also that all formulae here and in Appendix ( p 2[ ) hold 
for discrete space, under the same assumptions (defining Laplace transforms as B ± (u,u // ,p) — Ylt=i B ± (u, u" , x)), 
and one could also check whether discrete landscape could also be solved by these methods. 

Note finally that although fixed point (i.e. continuum) models can alternatively be solved by the constrained path 
integral method, the RSRG is in principle more powerful. Here for instance we only need assumptions (i) to (iii) to 
hold only asymptotically at large T. In principle the full convergence to the fixed point could be also studied, as in 
any RG method, and universal results established for a broader set of discrete models within the basin of attraction 
of the fixed points found here. 

In the following we restrict to the case of Schrodinger measures corresponding to the choice ( |67| ) and use the explicit 
solutions derived above to study the case of the toy model ([!]) in details and the case of stationary landscapes j2l| ) 
more briefly. 

IV. STATISTICS OF EXTREMA IN THE TOY MODEL 

In this section, we apply the general construction of the probability measure of renormalized landscape to the 
particular case of the toy model (Q). We start by giving the expressions for the blocks before constructing the full 
measure. 
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A. Blocks of the toy model in terms of Airy functions 



We have seen in ( |38| ) how the functional integrals for the process dll/dx = fix + rj(x) could be recast up to boundary 
terms into functional integrals of the Schrodinger-type. The Fokker-Planck blocks thus read (38) 



E± (u , uuxx- xo) = e Hu ,*o)-h(u uXl )E± (uQj XQt u ^ Xi) 
in terms of the Schrodinger blocks B E given in Laplace (j52|[36|) by 

1 



B T (u ,U,p) = B£(u,u ,p) 



k{u — r, u,p) 



exp( 



u-r 



dvdi In K(v, v + T, p)) 



E^(ub, u,p) = E T (u, itb,p) = exp( / di In K(u', u' + T,p)du') 



(69) 
(70) 



(71) 
(72) 



The function K defined in (|5^) now corresponds to the problem with the Schrodinger potential V[U] = (|4 
Two independent solutions of the Schrodinger equation are the Airy functions (Fl) 



4>i(u,p) — Ai(au + bp) 
(f>2(u,p) = Bi(au + bp) 



where 



1/3 



(73) 
(74) 



(75) 
(76) 



a is the characteristic inverse energy, whereas b represents the characteristic length scale of the toy model. The 
wronskian of the two solutions (|74|) reads 



and thus the function K reads(p3 



w{p) 



K(u, v,p) — — (Aifau + bp)Bi(av + bp) — Bi(au + bp)Ai(av + bp)) 
a 



(77) 



(78) 



Although there are no obvious simple expressions for the integrals involved in the blocks ([72]), in calculations it will 
be often useful to extract the leading behavior for u < v by rewriting 



K(u, v,p) = — Aitau + bp)Bi(av + bp) 1 
a 



and 



d\ lnK(u, v,p) = d u In Ai(au + bp) + d u In 1 — 



Bi(au + bp) Ai(av + bp) 
Ai{au + bp) Bi(av + bp) 

Bi(au + bp) Ai(av + bp) 



Ai(au + bp) Bi(av + bp) 



(79) 



(80) 



B. Normalization of the full Measure for the renormalized landscape 

As explained in Section (fl|), we consider left and right reflecting boundaries, i.e. we impose u(xl) = ul and 
u(xr) = ur with u(x2) = +oo and u(x^) — +oo. Then the normalization of the measure is ( |36| ) : 

Z[ul,xl;ur,xr) = / DUe 4 Ja= i Vda; ' = e *(*r-*l) (81) 

Ju(x L )=u L ^^(xr - X L ) 
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The decomposition into blocks at scale T reads 

+00 „ 2n— 1 



Z(u L ,x L ,u R ,x R ) — / J\ duidx i E^(u L ,x L ;u 1 ,x 1 )B r (u 1 ,x 1 ;u 2 ,x 2 ) 

™ — 1 J „• — 1 



n—1 i—1 

B+{u 2 n-2, %2n-2\ «2n-li ^2n-l)-^f («2n-lj ^2n-l5 ^H, (82) 
l-oo « 2n— 1 



e h(u L ,x L ) — h(u R ,x R ) 



^V^J n Alt -L 

^ / dUidXiE^ (ul,ui,xi - Xl)B^(ui,U2,X 2 - x%) . 

^ — 1 " ~- — 1 



n—1 i—1 

(w2n-2, W2n-l!^2n-l - X 2n -2)E^ (u 2n -l, UR, Xr - £271- 1) (83) 



i.e. the boundary terms involving the function h in the transformation from Fokkcr-Planck to Schrodinger ( p8[) cancel 
between two successive blocks, and thus only survive at the boundaries of the full system (xl, xr). 

C. Alternative expression of the full measure using barriers 

The invariance by a global shift in the energy u of the Fokker-Planck blocks 

Bp (tt ,ar ;wi,a;i) = Bp (0, x Q ; m - w ,^i) (84) 
B£(u ,x ;ui,xi) = B^(u - ui,x o ;0,xi) (85) 

translates into the following properties for the Schrodinger blocks ( |70| ) using the expression of the function h ( |39| ) 

Bf («o« "i.^i - a; o)= e ' l ( o . a; «)-' i («''^)+' l ("i^i)- ,l ("i-"«^i)i?-(0, It! - uo, an - x ) 

= e % xou °-% xlua B^ {0, Ul - u , Xl -x ) (86) 
B+(u , u u xi - x )= e Huo-uu*o)-Hu ,x )+h(u u x 1 )-h(a,x 1 )g+( Uo _ Ui xi _ x ^ 

= e* XoU1 -* XlUl B+(u - ui,0,xi - ar ) (87) 

and similarly for the edge blocks E. An alternative expression of ( |S3| ) containing only differences of energies is thus 
given by 

+oo r. 2n—l 

Z{u L ,x L ,u R ,x R ) = J2 / II duidxile^-^^-^E+iuL-ux^xr-XL)] 
n=l*' i=l 

Bp (0, u 2 - ui, x 2 - s 1 )e* X2 ( tt3 -" 1 ) J B+(u 2 - u 3 , 0, x 3 - ar 2 ) 
Bp(0, w 4 — ti 3 , x 4 — x 3 )e^ :E4 ^ 5 ~ u ^i3p (u 4 — us, 0, x 5 — ar 4 ) . . . 

Bf(0,u 2n ^ 2 - u 2n - 3 ,x 2n - 2 - x 2n - 3 )e^ 2;2 "- 2(tl2 "- 1 ~" 2 ' l - 3) B I !r(M2„-2 - u 2n -i, 0, x 2n -i - x 2n - 2 ) 
[e +f Kfl (» fl - U2n _ l) -^^Sp(0,un - u 2n _i,x fl - z 2 „-i)] 



where 



I f F - r 

B^(Q,F,p) = B+(F,0,p) = — — —— — exp( / <M> 4 In ^(u, u + I\p)) (89) 

E-{0,F, P ) = E+(F,0,p) =exp(+ / cfcdi lni<>, v + I\p)) (90) 

Jo 



D. Full measure using barriers for free boundary conditions 

The normalization ( |8l| ) concerns fixed boundary conditions {ul,Ur) and depends only on the difference (ur — Ul) 
since there is an invariance with respect to a global shift of all the energies u. For free boundary conditions, we need 
to sum over the end-potential ur with fixed ul, which of course yields a normalization to 1 
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/+oo 
du R Z(u L ,x L ;u R ,x R ) = 1 (91) 
-oo 

The decomposition into blocks at scale T now reads using the formulation in terms of barriers ( p8[ ) that we now 
denote by Fl = Ul — ui, F± = u 2 — u\, ... 

+00 r+oa r+00 f+oc 2 7t- 2 . 2 n— 1 

Z free (x L ;x R ) = l = J2 dFL dF R U dF * i I 

n —lJ0 JO JT i=1 J x L <x 1 <...x 2n -i<x R j—i 

[e^-^ F '-E+(F L ,0,x 1 - x L )] 

Br(0, F u x 2 -x 1 )e^ Fl - F ^B+(F 2 ,0, x 3 -x 2 ) 

SfT(0, - a^eW^-^Bj (F 4s 0, a* - a*) • • ■ 

Sr (0, F 2 „_ 3 , ^2„-2 - a; 2 „_3)e^ a:2 "- 2(i;,2 "- 3 - F2 "- 2) B+(F 2 „_ 2 , 0, x 2 „_! - :r 2 „_ 2 ) 
[ e +^BF R -4-lE-(0, F R , x R - x 2n . L )] (92) 



E. Full measure using barriers on the infinite line 



Of particular interest is of course the limit of the infinite system xl — > — oo and x R — > oo. We define the following 
function that describes the right boundary in the limit x R — > oo 



r+°° 2 
gr(x)= lim / dF R [e + ^ FR -^E^{0,F R ,x R - x)\ 



Using @ 



ST (a;) 



+ OC 



dX 
— ( 

2tt 



-iXx 



lim 



+00 



^a^+iArcR / ^jp^+^x R F ' — J Q F di/di In K(u,u+T,i\)du 



(93) 



(94) 



In the limit x R — > oo, the integral over F is dominated by a saddle with F large. Using (|8C|) it is convenient to rewrite 
<9i ]nK(F, F + T, iX)= d F \nAi{aF + biX) - ^ r (F, A) (95) 



aAi(a(F + T) + biX) 



irAi(aF + MX) [Ai(aF + biX)Bi(a(F + T) + biX) - Bi(aF + biX)Ai(a(F + T) + biX)} 



A saddle-point analysis with the asymptotic behavior (F5) immediately gives 



+00 



dC,e Cv Ai{C, + b) ~ e3 y 

y—Kx> 



As a consequence, we get 



9r(x) 



dX 
2V 



— zAa 



1 



aAi(ibX) 



e J o 



Similarly, the other boundary xl — > — 00 yields 

f +00 



lim 

XL — > — OO 



dF L [e^ x i-i XLF ^E+(F L , 0,x - x L )} = g r {-x) 



(96) 



(97) 



(98) 



(99) 



As a consequence, the full measure using barriers on the infinite line reads ( p2[ ) 

+00 /.-(-oo 2n— 2 „ 2n— 1 

Z /ree (-oo; +00) = 1 = ^2 / II ^ / IT dxjg T {-xi) 



Br (0, Fj.a^ - ii)e W^i-^B+^a, 0, z 3 - x 2 ) ■ • • 

Br (0,F 2n _ 3 ,a:2„-2 - a; 2 „-3)e*" 2 "- 2(F2 "- 3 ^ F2 "- 2) B+(F 2 „_ 2 , 0, x 2n -i - x 2n - 2 ) 
9r{x 2n -i) (100) 
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V. BOLTZMANN EQUILIBRIUM AT LOW TEMPERATURES OF THE TOY MODEL 



In this section, we use the explicit measure of the renormalized landscapes of the toy model derived in the previous 
section to compute various observables that characterize the equilibrium of the toy model ([!]) on the full line ] — 00, +00 [ 
at low temperature. For T = 0, we need to compute the distribution over samples of the position of absolute minimum 
energy on the full line. For small T, we need to compute the probabilities of configurations with nearly degenerate 
minima which govern the thermal fluctuations. 



A. Distribution of the position at zero temperature 

To obtain the distribution of the position of absolute minimum energy we need the final state of the renormalized 
landscape, corresponding to the r — > 00 limit, where all internal barriers have been eliminated. The decomposition 
( |100| ) then reduces to the term n = 1 that contains only the boundary blocks 

/+00 
dxg QO (-x)g QO (x) (101) 
-OO 

where (||) 

r+°° d\ e~ iXx 

&»(*)=/ ^-TT^ ( 102 ) 
J_ 00 2n aAi(biX) 

since the function (]9q ) vanishes for T — > 00. The distribution of the position x of the minimum on ] — oo,+oo[ 
reads 

P]-oo,+oo[{x) = g 00 {x)g 00 {-x) (103) 



in agreement with p7| , p8[ . The normalization is easily checked with the relation b = I /a 2 and the properties (F3) 
(0) 

+ °° dxP(x) = — - = h*^]*=+i°° = 1 (104) 

Let us now describe the asymptotic behaviors of this distribution. The Airy function Ai[z) has an infinity of poles 
which all are on the negative real axis : the largest one is e\ = —2.33810741.. where Ai'{— \e\\) = 0.70121082, the 
second one is e2 = —4.08794944.. where Ai'(— [eg]) = —0.80311137 and so on. For x < 0, we may deform the contour 
in the complex plane to express the function g as a series over the zeroes (— |e s |) of the Airy function 

+00 

g oo (x<0) = aJ2 TT, i (105) 
Thus g(x) decays exponentially for x — > —00. For x > 0, we may shift the contour integration A — > A — ic 

/+00 ^ — i\x — cx 

.„ ^ aM{UX + bc) (1 ° 6) 
Assuming c large, we replace the function Ai by its asymptotic behavior ( [F5| ) 

goo ( x ) = ' d\{b{c + i\)) 1 ' 4 e- x{ - c+ ^+i^ c+ ^ V2 (107) 



A saddle-point analysis leads to the choice c = and to the asymptotic behavior 



goo (x) ~ 2a 3 xe-^ a2x ^ (108) 

x — >+oo 



Putting (105) and (108) together, the asymptotic behavior of the distribution ( |l03| ) is given by 

1 _ 1 r„2ui-,3_ 21 



Pi-oo^x) c 2a 4 N , I'D'- l-ll-l (109) 
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B. Moments of the position at zero temperature 



Here we obtain the explicit analytical expressions for the disorder averaged thermal moments moments < x p > of 
the position in the limit of zero temperature, where they are determined by the absolute minimum. They can be 
derived from the Fourier transform 



/+oo <H 
dxe lsx P{x) = I 
-OO J — t 



- (110) 

2Tra 2 Ai(ib\)Ai(ib{\ + s)) V ' 



as 



+ °° du 1 

where C p are numerical constants given by 



<^r> = b P I ™_±_flj^ f ) = m - 2 p/ 3 C p (111) 



Cp = 2^ -J- flP ■ (-i-,) (112) 



For odd p, an integration by parts shows that C p vanishes, i.e. the odd moments vanish as it should from the symmetry 
(x, —x). For even p, we may also write 

We give the lowest moments: 

4/ , / ,+0 ° dfiii'H 2 2 4 / 3 /• +00 du -ra , . 

C 2 = -2 4 / 3 / ; ' = / — — = 1.05423856519.. 114 

2vr Ai(iu) 4 3 2tt^(w) 2 v ' 

Ci = 2 m = 3.15451.. (115) 

These exact results are in reasonable agreement with numerical simulations of Ref. [Q which found C2 = 1.05381 ± 
.01 and 3.15441±.04. It should also be compared with the result of approximate replica methods which give respectively 
C 2 = 3/(2Vtt) 2/3 = 1.29038 (RSB gaussian variational method), C 2 = 1.00181 (vector breaking) Q. Note also the 
exact result for the ratio R = C 4 /C| = 2.83827.. in good agreement with the numerical simulation result R = 2.84048, 
while the replica approximations give respectively R = 3 (GVM) and R = 2.995 (VB). 



C. Joint distribution of position and energy at zero temperature 

To obtain information on the minimal energy on the full line, we need to fix the absolute value of the energy at a 
given point for all samples. The usual natural choice is of course 

U{x = 0) = (116) 

Then in the renormalization procedure, we need to keep information on this value, that doesn't a priori appear 
in the list of extrema. As a consequence, to obtain the joint distribution P\-oo,+tx[( x 0i l *o) of the minimal energy 
U(xq) — —uo (with uq € [0, +oo[ and of its position xo, we need to consider the renormalized landscapes on the two 
independent half-lines ] — 00, 0] and [0, +00 [ which have the same statistical properties. More precisely, we have for 
the case xo > 

/•no r+oc 

P]-oo,+oo[{xo > 0,U ) = P[ 0t00 [(x ,Uo) dm dXiP[ 0iOO [(xi,Ul) (117) 

Jo Jo 

where Pm tO o[( x 0i ^o) is the joint probability that (xq,uq) characterizes the absolute minimum on the half-line, and 
where the integrals represents the probability that the minimum on the other half line is less deep than uq. The 



decomposition (92) for the case xl — and xr — > 00 becomes in the limit r — > 00 using (|93|) 
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Z /ree (0;+co) = 1 = / dF t 



+ 00 



+ OO 



dxEX 00 {F L ,Q,x)g 00 (x) 



where has been defined in ( 102 ) and where the edge block E + has for Laplace transform in the r — ► 00 

Ai(aF + bp) 



Ai(bp) 



As a consequence, the joint distribution for the half line reads 



P[o,oo[(x,u)= E$ oo (u,0,x)g oo (x) 



+OC 



dAi e lXlx Ai(au + biX^ 



2tt 



Ai(biXi 



+00 



dX? e 



-i\2X 



2-k aAi(biX2) 



in agreement with p7fl . In particular, for u — we have 

^[0,00 [{z,0)= S(x)g oo (0) 



(118) 



(119) 



(120) 



(121) 



as it should in terms of the function g^ ( jl02| ). 

For the distribution of the energy uo alone, the relation (117) becomes after taking into account the two sides x > 
and x < 



P]-oo,+oo[(uq) = 2P[o,oo[(wo) / duP[ 0t oo[{u) 



duo 



duP [0jOo[ (u) 



where 



+ 00 



P[o,oo[(u)= I dxP [0tOo[ (x,u) = lim 



+io ° dz 2 1 



2in aAi(bz%) 



dz\ Ai(au + bz\) 1 
2i7r Ai(bz\) p + Z2 — z\ 



(122) 



(123) 



is the distribution of the energy alone on the half- line [0, +00 [. For u > 0, we may close the contour in z\ on the right 
at infinity to obtain via Cauchy theorem 



-P[0,oo[(") : 



dzi Ai(au + bz 2 ) 
2iir aAi' 2 (bz2) 



+00 



dX Ai(au + biX) 
2^ aAi 2 (biX) 



(124) 



and this expression is also valid for u — (121). 
The expansion near u = reads 



P[n,oo[{u, 



+ OC 



dX Ai(biX) + auAi'(biX) + ^-Ai"(biX) 



2tt 



aAi 2 (biX) 



5oo(0)-y^(0) + O( U 3 ) 



which leads to the following expansion for the full line 



^-oo,+oo[(uo) ~ 2^05^(0) + O(u ) 



(125) 



(126) 



To obtain the asymptotic behavior for large u, we may shift the contour integration A — * A — ic, and assuming c 
large, we get via the saddle-point method with the choice c = |^ (F5) 



-P[0,oo[(") : 



+00 



dX Ai(au + biX + be) 



1 



27r aAi 2 (biX + bc) c->oo a 0r J_ 00 (au + biX + be) 1 / 4 



+OO 



dX- 



(biX + be) 1 ' 2 |ffra + 6 e )3/2 2 (oM+6iA+6e) 3/2 



~ —a(au) ' e 

u — >oo 3 



l /2 -4(^) 3 



(127) 
(128) 



D. Samples presenting nearly degenerate minima and thermal fluctuations at low temperature 

As explained in the introduction, there are remarquable identities Cq] for the disorder averages of the thermal 
cumulants of the position (^). It was moreover shown in |l8j that the result (J5J) for the second cumulant at low 
temperature contains a finite contribution from the rare configurations presenting two nearly degenerate minima. In 
this section, we compute the statistical properties of such configurations and show that they actually give the full 
exact result (S). 
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1. Configurations with two nearly degenerate minima 



It is shown in Appendix ( |d| ) that the probability to have two nearly degenerate minima separated by an energy 
difference AE — e — > and situated at x\ and x 2 reads 



T>(e,x 1 ,x 2 ) = egoc(-xi)d(x 2 - x 1 )g OQ (x 2 ) + 0(e 2 ) 



where the function g^ has been defined in (102) and where 



2tt Ai(ibX) 

The total probability density to have two degenerate minima separated by a distance y > thus reads 

+°° dX e~ iX y 



D(y)-- 



+0 °^ 1 lim/ P(e ' a;i ' a;i+y) 

— OO 



d(y) 



dxig{-x 1 )g(xi + y) = bd{y) 



2ir Ai 2 {ibX) 



In Laplace with respect to y we have 



D(p)= 



+ 00 



dye- py D(y) = a 



+io ° dz 2 1 



2iri Ai 2 (z 2 ) 



+ 200 



1 



dz\ Ai'(zi) 
2ni Ai(zi) bp + z 2 — z\ 



In particular, the second moment is given by 



dyy 2 D{y)= lim 



lim 



aV 



dp 2 



+lo ° dz 2 1 



+lo ° dz 2 



2iri Ai 2 (z 2 ) 
1 



2b 2 



2-iri Ai 2 (z 2 ) \ Ai 



too 
II 



dz\ Ai'(zi) 
2iri Ai(zi) (bp + z 2 — ziY 



(129) 
(130) 

(131) 

(132) 

(133) 
(134) 
(135) 



where we have used Cauchy theorem. This last integral may be computed via integration by part, with the help of 
the following properties of Airy function : the function f(z) — ^ife) satisfies f'(z) = z — f 2 (z) and the function 
= sati sfies g'{z) = -2f(z)g(z) : 



+oc ab 2 
dyy 2 D(y)= — 



+lo ° dz 2 1 



ab 2 



2iri Ai 2 (z 2 ) 



Bi(z) 
2iAi(z) 



z=-\-ioo 



Z — — IOQ 



2 



1 



(136) 



with (WB 



2. Contribution to the second thermal cumulant at low temperature 

The contribution of these configurations with two degenerate minima to the second thermal cumulant of the position 
(||) can be estimated at order T in temperature as follows. The two minima have for respective Boltzmann weights 
p = and (1 — p) — ]-f— -37- The variable (x— < x >) is (1 — p)(xi — x 2 ) with probability p and p(x 2 — xi) with 

probability 1 — p, so we get after averaging over e with a flat measure as in the study of the rare events in 



< (x- < x >) 2 >=p{l-p)(x\ -x 2 ) 2 (137) 

"+oo /-+00 /-+oo g — e/T 

dxi / dx 2 / dtV'(e = 0, xi, x 2 ) -e/t\ 2 ( ,T2 ~ ( 138 ) 

+ OG 

T I dyD{y)y 2 (139) 
o 



Using (13£), we exactly find the exact result (||). This shows that the thermal fluctuations at low temperature are 



entirely due the presence of metastable states in rare disordered samples. In particular, the susceptibility 
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^(<x 2 >-<x > 2 ) (140) 



has for average a finite value at zero-temperature 



X = ~ (141) 



but only the samples with two nearly degenerate minima contribute to this zero-temperature value. The typical 
samples with only one minimum have a susceptibility that vanishes at zero temperature. This type of behavior was 
also found numerically in |25| and analyzed in p0(] for the random directed polymer model. 



3. Contribution to the thermal moments at low temperature 
The small-T behavior of the even moments of (x— < x >) can be similarly computed 



< (x- < x >) 2n >= - / y 2n D{y) + OiT 2 ) (142) 
n Jo 



The moments of D(y) read, with the same method as in (135) 



^DM-^-DV/ 5SH y M-jpputi- (143) 



where vrik are numbers given by 



mi = ( _ 1 )^-«/3/ + ^_i_(^) < " (2) (144) 



4- Discussion 



The comparison between (142) and (Q) shows that there are a lot of cancelations in the disorder averages of the 



cumulants. This phenomenon was already emphasized in |19| | for the random directed polymer, and will be studied 
in more details in ]39[ | for the toy model as well as in other one-dimensional random systems which do not present the 
statistical tilt symmetry leading to the identities (^). 

Another question that needs further explanations and that will be studied in |39| is why the result (|5|) which can be 
understood at order T by considering only the configurations with two degenerate minima remains true for arbitrary 
temperature. Indeed, in the light of the above discussion, it would be natural to expect corrections at order T 2 coming 



from configurations with three degenerate minima and described by the distribution function analogous to (129) 



T(ei,e 2 , xi, x 2 ,x 3 ) = eie2goo(-xi)d(x 2 - ^1)0^3 - x 2 )g oc (x 3 ) (145) 

and so on for configurations with n nearly degenerate minima. One also expect corrections from the thermal width 
around minima. The study of the roles played by these various effects in different systems is postponed to a future 
work il. 



VI. NON-EQUILIBRIUM DYNAMICS OF THE TOY MODEL 

In this section, we study the non-equilibrium dynamics of the toy model by using the full measure of the renormalized 
landscape of the toy model for an arbitrary scale T. We also briefly summarize some results for the RFIM in a field 
gradient, already discussed in the Introduction. 

As in the Sinai model the long-time dynamics is asymptotically entirely determined by the renormalized land- 
scape at scale 

r = Tln< (146) 
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in the sense that a particle starting at xq at time t — will be at time t around the minimum of the renormalized 
valley at scale T that contains the initial position xq. As in the Sinai model this method allows to compute 
asymptotically exact results for scaling variables in the limit t — > oo corresponding to large barriers, and to large 
distances. Here in the toy model, there is a typical scale of energy given by (1/a) and a typical length b = 1/a 2 fl76|). 
As a consequence, we need to define the following scaling variables for the barrier 

j = aT=(^) 1 / 3 T\nt (147) 

and for the length 

* = f = (f ) 2/3 * (148) 

The RSRG method will thus give asymptotic exact results for the dynamics in the double limit t — ► oo (to have 
large barriers T) and /i — * (to have large typical scales for energy and length in the toy landscape) with the scaling 
variables 7 and i fixed. For large 7, as studied in previous section, the system will reach equilibrium. For small 7, i.e. 
when the barrier T = Tlniis much smaller than the characteristic energy scale (1/a) of the toy model, the dynamics 
should look like the Sinai model, where the typical length associated to the dynamics is T 2 Q so that the scaling 
variable for distance of diffusion during the time t reads in this regime 



'A-Yk =^ = ^ I 1-1!) I 



We will first consider a global quantity of the whole system : each sample is characterized by an equilibration time 
t eq corresponding to the largest barrier T eq in the system via 

T eq =T\nt eq (150) 

We will compute the distribution over samples of this equilibration time and then characterize the samples which 
have an anomalously large equilibration time and which thus dominate the long time dynamics. We will then study 
local quantities, such as the probability density of minima of the renormalized landscape, and the diffusion front, i.e. 
the generalization of the Kesten distribution of the Sinai model to the case of the toy model. 

As discussed in the Introduction, the equilibrium state for the RFIM in a field gradient for large J is obtained by 
stopping the RSRG at scale T = T eq o = Thit eq o — 2 J at which equilibrium is attained and putting a B domain wall 
in each minimum and a A domain in each maximum of the renormalized landscape at scale T = 2 J. The equilibrium 
state thus consists in a B domain plus n additional pairs of domains A — B, the probability of n with respect to 
sample being thus N^i^T^ ■ The distribution of the total number of domains at equilibrium is thus a scaling function 
of "feqo = 0,2 J. Below we obtain the behavior at large "f e qO of these probabilities for n = 0, 1 as well as an estimate of 
the total number of domains at small 7 e qo- Similarly, inserting instead T = Tint, one also finds these probabilities 
for the out of equilibrium system at time t starting from a random (i.e. high temperature) configuration at time zero. 
The total magnetization reads 

N A N B 
a—l a— 1 

its probability measure at time t can thus be formally expressed using the full measure for the renormalized landscape 
at r = Tint, and the equilibrium expression obtained using T = 2 J. In the limit of large jeqo, its probability 



distribution coincides exactly with the one of the absolute minimum of the toy model given by formula (103) setting 
M = 2x. 

For non random initial configurations prepared with a single domain B 1 the equilibration time of this domain is t eq 
computed below (provided Tlnt eq o > T\nt eq , otherwise additional pairs are created at time t eq o and take positions 
in the additional deep extrema). If the single domain is A, then it is expelled from the system when the extrema 
corresponding to the highest barrier is decimated, again the expulsion time is the same T\nt exp = T\nt eq . Note, 
interestingly, that if two almost degenerate minima (within energy difference T) exist, the domain wall A will be 
expelled with finite probabilities towards the two boundaries, and the system will be in a macroscopic mixture of two 
reversed states (either all spins (+) or all (-)) between times t exp and t eq o (if we assume that hL <^ J). 
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A. Distribution over samples of the equilibration time 



We first study the distribution over samples of the equilibration time t eq at which equilibrium is reached. The 
probability P(t eq < t) that the system has already reached equilibrium at time t is directly given by the probability 
that there remains only one valley at renormalization scale r = Tint in the decomposition (100) with (|9q) 



r+oo r+oo ^ 

P(t eq < t)= I dxgv(-x)g T {x) -- 



1 



2tt a 2 Ai 2 (ibX) 



e -2j+°°dFi>r(F,\) 



(152) 



where tpr has been defined in ( |96| ) 

e it as . 

V(t eq <t) = $( 7 ) 



It is convenient to rewrite it as a scaling function of 7 defined in (147) 

+00 



with 



V> 7 (/,A) = 



dX 1 p -2j +°°rf/^(/,A) 

2irAi 2 {iX) 
Ai(f + j + iX) 



nAi(f + iX) [Ai(f + iX)Bi(f + 7 + iX) - Bi(f + iX)Ai(f + 7 + iX)] 



(153) 



(154) 



Equivalently, this result describes the probability distribution P ma x(F) over samples of the largest barrier T 
T In t eq existing in the system 



(155) 



1. Probability of large equilibration times 



For large 7, the dominant behavior of the function i/> 7 (154) is found by using ( |F5| ) 

Ai(j + f + iX) : 



which leads to 



_ e -|(7+/+^A)) 3 / 2 

Si (7 + / + iX) 7-00 2 

ae -|(7+/+jA)) 3 / 2 



2irAi 2 (f + iX) 
1 

7^00 27rAi 2 (/ + iA) 6 
The change of variable / = gives 

'•+ 00 dw; 1 





-iA) \] 




HA) J J 



| 7 3/2_ 27 l/2 (/ + lA) + ^_l_^ 



7-~y V7 2^2^A+-^) 



e _| 7 V2_ 24A7 l/2_ 2u) 



1 + 



,1/2 



1 



7 — >oo 47r v / 7^4i 2 (zA) 



| 7 3 / 2 -2iA 7 1 / 2 



1 + 



yl/2 



The leading behavior of ( [15 3D is thus given by 



$(7) 



+OO 



dX 1 



7— >oo j /_ 00 27rAi 2 (iA) 

_ 4 
e 3 



1 



1 



7^00 2^^ J_ OQ 2tt Ai 4 (iA) 



2 7 r % /7^ 2 (iA) 
1 + 



e _| 7 3/2„ 2jA7 l/2 



1 + 



yl/2 



yl/2 



Shifting the contour integration A — » A — ic and assuming c large, we get using (F5) 



(156) 

(157) 
(158) 

(159) 
(160) 

(161) 
(162) 
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dX e ~ lX2 ^ f +oc dX e -Vr(«=+<A) 



2tt Ai 4 (iX) 



2ir Ai 4 (iX + c) c^oo 



87T 



+ 00 



dX(c + iX)e- 2 ^ c+ ^ + ^ c+lX ~> 3/2 



A saddle-point analysis leads to the choice c = ? and to the asymptotic behavior 



(163) 



.oo 27r Ai 4 (iA) 7^+00 



^V2^7 5 / 4 e"^ 3/2 



(164) 



The probability to have only one valley at time t finally behaves for large 7 as (162) 



V{t eq <t) = $(7) 



7 3/4 e -|7 



3^3/2 



(165) 



Equivalently, this result describes the asy mpt otic tail of probability distribution P maa; (r) over samples of the largest 
barrier V = T\nt eq existing in the system (155) 



5/4 -|(or) 



3/2 



(166) 



2. Small equilibration times 



For small 7, the leading behavior of (153) can be found with the change of variable / = ^ using (F5) 



+00 



2sds 



Ai^+iX + j) 



T Jo 



+00 



dss* 



1 n Ai(^ + iX) \ Ai($ + iX)Bi(^ + iX + 7) - Bi{$ + i\)Ai($ + iX + 7) 
C(3) 



sinh s 7 3 



(167) 



(168) 



where C( n ) is the Riemann zeta function. This leads to the following essential singularity for the short equilibration 
times 



V(t eq < t) = $(7) ^ e ^ 
7— >o 



(169) 



or equivalently for probability distribution P ma x(^) of the largest barrier T existing in the system (|155| 



r m „ T 1 = — <P al ~ -r-e <° r ) 

v ; dr v ; r^o a r 4 



(170) 



B. Rare configurations with anomalously large equilibration times 



We are now interested into the statis tical properties of the configurations that present an anomalously large equi- 
libration time, as described by the tail ( 165 ) since these configurations are important f or th e long time dynamics. In 
the regime 7 > 1 under consideration, the leading correction to the one valley event (16E) that we have computed 
in the previous section corresponds to the two valleys configurations, because the configurations with three or more 
valleys are even more unlikely. These configurations correspond to the last non-equilibrium dynamics in the system 
before it finally equilibrates by going over the largest barrier T eq . We can already deduce from the previous section 
the total probability of the two valleys configurations (165) 



N, 



3/4 e -f(ar) 3 



(171) 



To characterize the statistical properties of these two valley events, we have the probability measure (100) 
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N r 4) {x 1 ,F 1 ,x 2 ,F 2 , x 3 ) = g r {- Xl )Br(0, F 1 ,x 2 - x 1 )e^ Fl - F ' ) B+(F 2 , 0, x 3 - x 2 )g r (x 3 ) 



(172) 



where x\ < x 2 < x 3 , X\ and x 3 are the two local minima separated by the barriers F% > T and F 2 > T from the local 
maximum x 2 . For large T, the leading behavior of the blocks B (BOh can be found using (79) 



K(u,v,p) ~ — Ai(au + bp)Bi(av + bp) 

v>>u a 



which yields 



Br(0,F,p) = B+(F,0,p). 



r~*oo TrAi(bp)Bi(aF + bp) 
At the same level on approximation, t/p is given by g^ ( |102| ) , so that we finally get at leading order 



Np(x u F l ,x 2 ,F 2 ,x 3 ) 



e^xaCFx-ft) /-+OC ,+oo ,+oo ^ ,+oo 



2tt 



2tt 



-oo — oo 

:J i(Ai + A2)2;i+i(A3-A2)2;2-j(A2+A3)a:3 



2tt 



2tt 



Ai(WA 1 ) J 4i(6A 2 )-Bi(aJ , 1 + bX 2 )Ai(bX 3 )Bi(aF 2 + bX 3 )Ai(biX 4 ) 



(173) 
(174) 

(175) 
(176) 



It is shown in Appendix (^) that after integration over the barriers and the position x 2 of the local maximum, the 
measure for the two minima (x i,x 3 ) has the simple structure 



M^\ Xl ,x 3 y. 



/ dx 2 / dF 1 / dF 2 N^\x u F u x 2 ,F 2 ,x 3 ) 
Jx! JT jt 

r5oo 2 ^( a a r)V2 e ' |(ar)3/2 g°°^ :Cl ^ ) ( a:i ~ x * + nar) 1/2 ) 9oo (x 3 ) 



with 



4-oo 



dX 



2tt a 2 Ai 2 {biX) 



{goo *goo)(x) 



(177) 
(178) 

(179) 



The total normalization (171) is recovered by integration over (xi,x 3 ) with the use of (164). 

Since all the functions 3oo(y) are centered around y — 0, the two-minima measure M^\x\, x 3 ) is concentrated 
around (xi,x 3 ) — (~xr,xr) with 



X T 



2a 2 



(180) 



More precisely, a saddle-point analysis leads to the approximation by a product of two Gaussian distributions of 
broadness 



V3 



as 



2a 2 (aT) 1 / 4 



(4) e e 2 "r 



M r 4) ( Xl ,x 3 ) ~ = JV£ 

r V ^r^co r 2v/2^ 2^2^ 



(181) 



(182) 



concentrating all the normalization found above (171). 

The two valleys configurations with an anomalously large barrier ar > 1 thus present minima which are anomalously 
far from the origin when compared to the typical scale of minima (|76|) 



(183) 



c y ! = > b 



which is necessary with the Brownian scaling to have a barrier of order T between them. 
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C. Probability density of minima in the renormalized landscape 



The probability Vr(x) that the point x is the bottom of a renormalized valley at scale Y can be decomposed into 

Vr(x) = V T L (x)V r R {x) (184) 

where V R {x) represents the constraints on the landscape on the half line [x, +oo[ : the potential starting from U(x) 
has to reach U(x) + Y at an arbitrary position y without any return to U(x), which reads in terms of the Fokker-Planck 
block with the use of (USt!) 



V R (x)= 



+ 00 



dyB r (0,x;T,y) 
+OQ dX 



+OO 



dye h ^- h ^B-(0,T;y - x) 



-iXx 



, 2tt tt {Ai{ibX)Bi(aY + ibX) - Bi(ibX)Ai(aY + ibX)) 
so that finally in terms of the scaling variables 7 ( |147D and x ( ]148| ) 

yR( )— Si 3 s / - 

r W- - e3 J ^ 2^ (Ai(iX)Bi(j + iX) - Bi(iX)Ai(j + iX)) 

Similarly, for the left half line ] — 00, x], we have by symmetry 

V r L (x)=V R (-x) 



dye lXy+ ^ r -^y 3 



dye 



(185) 
(186) 

(187) 



(188) 



In the limit 7 — > 00, we have using ( 173| ,F5) and by evaluating the integral over y by the saddle-point method 

ho ° dX e~ lXi r+ °° 



7 R (x) ~ -e^ 3 

7— >oo 7T 



2tt Ai(iX)Bi(~i + iX) 



dye 1 



+ °° dX e- lXi 



2tt Ai(iX) 



is 3 / \ 

e 3 5oo(x) 



leading to 



V^x) = g c(~x)g 00 (x) 



as it should to recover the probability distribution of the absolute minimum of the full landscape ( 103 ). 
In the other limit 7 — > 0, we have using (F2) 



V r R (x) ~ ae* £ 
7^0 



+00 



dX 



2tt 



(; • -'A • •••) 



+00 



dye 1 ^-^ 3 (1 + 7J/+ ...) 



7^0 7 



(1 + 7X) + 0(7) 



and thus, for the points x that are in the region where x remains finite, the probability (184) is given by 

V r (x) ^ n 4 = 4 
7^0 7 Z l z 



(189) 
(190) 

(191) 

(192) 
(193) 

(194) 



i.e. it represents in this regime as it should the uniform density 1/r 2 of minima of the Sinai model Q. 

Another quantity of interest is probability that x is the nearest minimum x^ in from the boundary x R = +00 in 
the renormalized landscape at scale Y 



Wt?(x) = V^(x)V^(x, +00, +00) 



(195) 



where the constraint on the left half line is the same as before ( 184. 188] ) and where the constraint on the right half- line 
reads (p|) 



+ OC 



V R (x, 00, 00) = lim / dF R E T (0, x; F R , x R ) = e^ x g r (x) 



(196) 
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or in terms of of the rescaled variables using (|9£ , 154) 

" +0 ° dX 
, 2? 



y r ^(x, oo, oo) = ae» 5 



— iAa 



(197) 



Of course, in the regime 7 — > 00, the distribution (195) of x^ in is the same as the distribution of the absolute 
minimum (191). In the other limit 7 — > 0, we expect that x^ in is much bigger than the scale b (f76j) of the absolute 
minimum. In this regime, its position will not be governed by the constraint on right landscape Vfrt (x, 00, oo), which 
is easy to satisfy in this region, but by the constraint on the left landscape, which behaves in this regime as ( |193| ) 

V r L (x) = V r R (-x) ~£(l_ 7 £)+ ( 7 ) 
7^0 7 



(198) 



We thus expect that x^in sca l es as l/7> i- e 



' 7^0 7 



(199) 



As a consequence, the total number n of minima is the renormalized landscape at scale T can be roughly estimated 
in this regime by dividing the effective region (199) by the typical distance T 2 (194) between two minima 



-.typ 



7^0 7 



(200) 



This estimation is coherent with the behavior found in (169). Of course, an exact calculation of the distribution of 
the number n of minima of the renormalized landscape could in principle be done from (100) 



Z free 1 



+00 
n=l 



(2n+2) 

r 



(201) 



where N^ n+2 represents the total normalization of the renormalized landscapes containing exactly n local minima. 
We have considered previously ( 152 ) corresponding to a single minimum, and ( 171 ) corre sponding to two 
minima, but unfortunately, the multiple integrals over the barriers F and positions x needed in ( |l00|) become rapidly 
intractable for arbitrary n and arbitrary 7. 



D. Diffusion front 



As explained at the beginning of the Section [v|, the long-time dynamics is determined by the properties of the 
renormalized landscape T = Tin t, an d the re sult s obtained will be valid in the double limit t — > 00 and /1 — > in 



terms of the scaling variables 7 (147) and x (|14S| ). As in the Sinai model ||, the average over the disorder of the 
diffusion front P(x, t\xo, 0) is given by the probability that x is the minimum of a renormalized valley at scale T that 
contains the initial position xq. 

Le t us first consider the case x < xq. The constraint for the left half-line ] — 00, x] is simply given by Vt (x) = Vr{— x) 



(184,188). The constraint for the landscape on the right half line is now that x should be the starting point of a 



renormalized bond of arbitrary barrier F > T and of length I satisfying x + I > xq. We thus have 

/ r+00 r+00 \ 

6(x -x)P(x,t\x ,0) = Vr{x)[V-r R {x,+oo 1 +oo)+ dF dlV R (x,F,l)) (202) 

V Jr Jx ~x J 

where Vp represents the constraints on the landscape o n th e half line [x, +00 [ : the point x should either be the first 



minimum on the right of the renormalized landscape ( 196 ) or it should be the starting point of a finite ascending 



renormalized bond of barrier F and length I, followed by an arbitrary descending bond, i.e. in terms of the Fokker- 
Planck block ©H) 

r+oc 

Vj?-(x,F,l)=Br(0,x;F,x + l) / dyB+(F, x + I; F - T, y) (203) 

Jx+l 
r+00 

= Br(0,x;F,x + l) dyB+(T,x + l;0,y) (204) 

Jx+l 
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Using the explicit expressions (^,^,^8|), we have 

y«( : r,F,0=4 e " x3+f(F " r)(a;+ ' ) 



7T 

+00 



dXi 



JMI + J CtVa{ A,{av + b,\ 1 )B,( a (v + r) + bi>. 1 )-Bi(av-l-bi\ 1 )Ai{a{v + r) + b,X 1 ) ) 



2tt (Ai(a(F — T) + biXi)Bi(aF + MXi) - Bi(a(F -T) + bi\i)Ai(aF + bi\ x )) 



+00 



dye " v 3 



+ °° dX 2 



,i\-2(y—x—l) 



, x+l 2tt (Ai(bp)Bi(aT + bi\ 2 ) - Bi(bp)Ai(aT + bi\ 2 )) 

so that in terms of the scaling variables 7 ([1.47D and x (|148|), the integral needed in (|202|) reads 



+30 



dF / dZ Vr(x,F,Z) 

00 z'+og 



00 IP r+00 



df 



bdl Vj?(x = bx,F = t,l = bl) 



df / dUeU-rW) 



a 1 z.3 

7T 

,X r. r/— r J-.Y A,'( 1 »+it 1 )B,(»,+ 7 +it 1 )-J3i'[» + iA 1 ) J 4.(i.+ 7 + i> 1 ) ' 

+ °° dAl e +J ° d ™( Ai( lu + iA 1 )B I (™ + T + iA 1 )-Bi ( ™ + iA 1 )>l I (™ + 7 + iA 1 ) , 



2tt (,«(/ - 7 + *Ai)B*(/ + *Ai) - Bi(f - 7 + *Ai)Ai(/ + iAi)) 

e iX 2 (y-x— I) 



- - dX 



dye 



x+l 



2tt {Ai{iX 2 )Bi{ 1 + iX 2 ) - Bi(iX 2 )Ai(j + iX 2 )) 



(205) 

(206) 
(207) 

(208) 
(209) 

(210) 
(211) 



In the limit 7 — > 00, the term (196) dominates and the diffusion front converges as it should towards the distribution 
of the absolute minimum for x and doesn't depend anymore on the initial position Xq 



P(x,t\xo,0) ~ g oc {-x)g 00 (~x) 

7 — >oo 



(212) 



In the other limit 7 — > 0, the term (196) is negligeable, and the typical length associated to the dynamics is T 2 
(|194|) . The diffusion front is expected to become in this regime a function of the scaling variable (|149|) 



(AX) 



Sinai 



Xq — X Xq — X 



r 2 



r 



(213) 



The contribution of finite barriers (211) indeed becomes with the appropriate changes of variables I = 7 /, Ai = -4-, 
I = 7(1 + V) 



+00 r+oo r+oo r+00 

dF dlV R (x,F,l) ~ -e* 5 / dTj / die 7 " (£+7 



riA 



iA;+ /J 1 7du 



Ai'(T J +i4)Bi( 7 »+ 1 +iA)-fli'(7» + iA)Ji(7»+7+iA) 

2_ 21 11 11 

Ai( 1 «+ii)B.( 7 «+7 + iA)- B, h»+.4-) Ai h» + 1 + '4) 



2vr 



(^(777 + t^)B*(7(l + r?) + i£) - Bi^rj + i^)Ai(7(l + 77 + 



-i(cE+7 2 3 



(214) 

(215) 
( 216 ) 

7^0 7 Jo 

where we have recognized the joint distribution of rescaled barrier 77 and rescaled length I of the Brownian landscape 



-00 r+oo 

di] / dlP*(r),l) 



p*{n,i) 



+00 



1 /2 

^ A (j^j e a;- 7) (a) 1 / 2 coth((a) 1 / 2 ) 
2tt S inh (iA) 1/2 

The diffusion front has thus for Laplace transform in this regime 



(217) 



"+00 r+00 r+00 r+00 1 _ „— pi" 1 1 

dioe-^-^Pfx.tlio.O) ^ / dye- pT v dl dr,P*{r,,l) = -=- (1 - coshr^p) (218) 

v ■y^°Jo Jo Jo p t 1 p 
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which is exactly the Laplace transform of the Kesten distribution for the Sinai model 
The diffusion front for xq < x can be obtained from the symmetry x — > — x 



r+cc r+oo 

6(x-x )P(x,t\x ,0)= [Vr(x,+oo,+oo)+ dF dl V L (x, F, I) ) V*{x) (219) 

V Jr Jx-xq J 

/ r+oo I-+0O \ 

= ( V^{-x, +00, +00) + / dF dl V R (-x, F, 1) ) Vf(x) (220) 



VII. RESULTS FOR STATIONARY LANDSCAPES 



In this section, we discuss the case of the stationary landscapes defined by (|T]) which belong to the class of solvable 



models with our method as discussed in Section (III A) 



A. Explicit expressions for the blocks 

We have seen in (|4l| ) how the functional integrals for the process dU jdx = —dW[U]/dU + r](x) could be recast into 
a functional integral of the Schrodinger-type up to boundary terms. The Fokker-Planck blocks read ( ^l|) 

^(uo.io.m,^) = e-& w M- w to»B£(uo,xo,ui,xi) (221) 
^(tto,io,m,n) = e-^ w ^ w ^E±(u ,x ,u llXl ) (222) 

in terms of the Schrodinger blocks B and E given in Laplace transform by 

1 r~ r 

B T (u , u,p)= B£(u,u ,p) = — rexp(/ dvd x hxK(v, v + T,p)) (223) 

(u,u b ,p)= (u b ,u,p) — exp( / dvd\ lnK(v, v + T,p)) (224) 

J u 

The function K (|53|) now corresponds to the problem with the Schrodinger potential V[U] = — ^W"[U] + jW'[U] 2 . 

For the special case p = 0, two independent solutions of the Schrodinger equation are then known for an arbitrary 
potential W[U] 

Mu) = e-? w(u) (225) 
<h(u) = e~ iw(u) C e w ^du' (226) 







of wronskian w = 1 and thus the function K ([53|) simply reads : 

K(u,v,p = 0) = e-^ w W+ w W I e w ^du' (227) 



The Schrodinger blocks of the measure read (224) 



Br(uo,u,p = 0) = B+(u,u ,p = 0)= * ' w * e J "° /" +r i^™W (228) 

Ju— r e u 

E+(u b ,u,p = 0)=E-(u,u b ,p = 0) = e 2 2 Ju ^^^W (229) 
so that finally the blocks for the full measure (E23) read in terms of the arbitrary potential W(u) 



29 



e W(uo) 



^ l»-»-/-")= p \ W{u , )du e !: ° /^--— (230 ) 

J u—T 



B+(u,u ,p = 0)= f: J eW - ul)dui e-^ ,.-.™-wcw (231) 

= 0) = ^W-^^-" 1 T^I^wt^ (232) 

(u 6 , u, p = 0) = e" /r 1 '^ du , (233) 

These expressions for the Laplace variable p — describe the properties of the extrema in the energy u of the landscape 
after integration over the positions variables x. 

B. Normalization of the full measure 

The normalization for left and right reflecting boundaries, corresponding to u{xl) — «l and u(xr) = ur with 
u(x2) — +oo and u(x^) = +oo, is given by ( |3^ ) 



W(u) _ ru 

■I 11 1 



Z(u L ,x L ;u R ,x R )= DU exp(- / dx(-(— - F[U}) 2 + -F'[U] j) (234) 

J u(xl)=ul Jxl 



and the decomposition into blocks at scale T reads ( |30| ) 

+oo 2n— 1 

Z(u L ,x L ,Un,x R ) = ^ / du i dx,^(uL,a;L;ui,a;i)S r :(ui,a;i;u2,X2) . . . (235) 

n=l J i=l 

■ ■ • -Bp ( u 2n-2, x 2ra _ 2 ; u 2 „-i,X2„-i)-B r 7(n 2 „-i, x 2 „_i; ur, x^) (236) 

+oo „ 271 — 1 

= e -i(w[u(x R )]-w[u(x L )])J2 / [] duidar^Cuz.ui.^-XijB^Cui.ua^a-a:!)... (237) 

n=l"' i=l 

Br(u 2 n-2,U 2 n-l,X2n-l - X 2 n-2 )Ef (u 2n - 1 , UR, XR - X 2 n-1 ) (238) 



C. Example : Brownian landscape with drift 

We rewrite the general expressions for the particular case W(u) — fu corresponding to the Schrodinger potential 

V(u) = f/4 = S 2 (239) 

where we have introduced the notation 5 = //2 to compare more straightforwardly with the papers @. The functions 
(j>i(u,p) and 4> 2 (u,p) introduced in (p3|) and their wronskian w(p) read 



0i (u) = e -"VP+^ (240) 

2 (tt) = e"V^ (241) 

to(p) = 2Vp+^2 (242) 

Thus the function if introduced in (153) reads 



This leads to the expressions 



lf(u,t>,;p) = pi smh{{v-u)^p + 8 2 ) (243) 
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B-( Ul ,u 2 ,p)= e-*(^-«0 ^ P + 62 exp(-( M2 - Ul -T)^/p~TP coMY^/p~P)) (244) 

sinh(ry?J + <5 2 ) 

= e -vs Vp + S 2 exp (_( U2 _ Ul _ r ) ((5 + ^+~p coth(rVpT^))) (245) 

smh(T y / p + S 2 ) 



B+(«i,«a,p)= e- 5 ^-" 1 ) — Vp + ^ 2 _ U2 _ r) coth(rv^T^)) 

sinn(r^/]3 + o 2 ) 

e T6 y/p + S 2 exp( _( Ml _ U2 _ r)(-5 + + <5 2 coth(rv5 + S 2 ))) (247) 

sinh(r N /p + (5 2 ) 



(246) 



£#"(ul, u)= exp(-(u £ - u)(-S + y/p + 5 2 cath(Ty/p + 6 2 ))) (248) 
Er(u,u R )= cxp(-(u R - u)(6 + y/p + S 2 coth(lVp + <$ 2 ))) (249) 

in agreement with |T^| . We thus recover simply and from a more general perspective the RSRG fixed-point for the 
pure Brownian landscape . 

D. Some further examples 

One can study cases where the potential W(u) (|2l]) is a simple power law W(u) — ru 1 . This corresponds to the 
Schrodinger potential ( [42[ ) 

v[u) = _n(7-i) u7 - 2 + i^! u 2 (7 -i) (250) 

which a priori contains two power-laws. Apart from the case 7=1 corresponding to the biased Brownian motion 
discussed above ( |239| ), there are two other simple cases 7 = 2 and 7 = 0, which we now briefly comment on. 

1. Landscape produced by a Brownian particle in a parabolic well 

The parabolic potential W(u) — ru 2 corresponds to the quadratic Schrodinger potential 

V{u) = -r + r 2 u 2 (251) 
The functions 4>i(u,p) and (j>2(u,p) introduced in fl5Sj ) and their wronskian w(p) read 

M^P) = e-^^uU{\ + |r U 2 ) (252) 

fa(u,P) = e-^ 2 uM(± + |-, |™ 2 ) (253) 

»(p) = rf^ljL) ( 254 ) 
where U (a, b, z) and M(a, b, z) are the confluent hypergeometric functions. In particular for p = 0, we have , § , z) — 



l/x/z and ( |226| ) yields 

&(u,p = 0) = e~^ =e -i^(«) (255) 
2 (u,P = 0) = e~i" 2 / dw'e™' 2 (256) 







u;(p = 0) = 1 (257) 
K(u, v,p = 0) = e -i {u2W) { du'e ru ' 2 (258) 



The above expressions allow to obtain the full joint distribution of wide excursions (extrema) of the trajectory U{x) 
of a Brownian particle in a quadratic well, but we will not explore these results further here since this goes beyond 
the scope of this paper. 
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2. Landscape produced by a Brownian particle in a logarithmic well 



We now consider the logarithmic potential W(u) = — rhxu for u €]0, +00 [ corresponding to the Schrodinger 
potential 



(TV 

The functions 4>i(u,p) and <f> 2 (u,p) introduced in (|5^) and their wronskian w(p) read 

<j>i(u,p) = yfuKi^{yJpu) 

<p2{u,p) = Vull^zi^/pu) 

w(p) — 1 



where K v {z) and I u (z) are the Bessel functions. For the special case p — 0, (|226|) yields for the function K 

1 



K(u,v,p = Q)=e-^ w ^+ w ^ / e w ("'W = 



1-r 



u r/2 v l-r/2 _ ^r/2^/2 



(259) 



(260) 
(261) 
(262) 



(263) 



E. An example with a long-ranged constraint : Positive Brownian landscape 

As an example of long-ranged constraint, we consider the process U{x) = |^(a;)| where V(x) is a pure Brownian, 
i.e. U(x) is constrained to remain positive by a reflecting wall at U — 0. For fixed boundary conditions (ul, xl) and 
(u Rl x R ), the total normalization reads 

Z(u l ,xl',u r ,xr) = e ^r-vl) +e 4 (*r~*l) (264) 



yjAir(x R - x L ) 

In Laplace transform with respect to the length (x R — Xl), this leads to 

r+00 1 . . 

Z(u L ,u R ;p) = / dx R e-^- x ^Z(u L ,x L ;u R ,x R ) = — (e-^ UR - u ^+e-^ UR - u ^) (265) 
Jx L 2y/P V ' 

We again consider the convention U{x~ R ) = U(x^) = +00 to define the renormalization procedure. The decompo- 
sition of the renormalized landscape at scale T can then be written as 

+00 . 

Z(u L ,u R ;p) = E^(u L ,u 1 ;p)B^(u 1 ,u 2 ;p)...E^(u 2n ^ 1 ,u R ;p) (266) 

n=l Ju i>° 

+E+ (u L , 0; p)G r (0, 0; p)Ef (0, u R ; p) (267) 
+E+(u L ,0;p)Gr(0,0;p)\Y / I B^{0,u 1 ;p)B+(u 1 ,u 2 ;p)...B+(u 2n ^ 1 ,0;p)\ G r (0,0;p)E^ (0,u R ;p) + ... (268) 

\n=l J «i>« / 

wher e the blocks {B,E) represents the 'bulk' blocks and are thus identical to the blocks for the pure Brownian motion 
( p9|) for (5 = 

B-(«i, ua; p)= B+{u 2 , Ul ;p) = 0((u 2 - u x ) - r) . -^^ -((^-"O-r^cotMryp) (26Q) 

smh(I y/p) 

Ef (u, u b ;p)= E+ (u b ,u; P ) = 8(u b - u ) e -K-«)v^oth(ryp) (270) 

and where the notation Gr(ui,u 2 ,p) represents the measure of the Brownian paths that remain in [0, T] between 
U\ and u 2 . Since U{x) = \V(x)\, this is equivalent to asking that the unconstr ained Br ownian V(x') remains in the 
interval ] — T, T[, and goes from itj to (+u 2 ) or to (—u 2 ). With the notations of (CI C14), its Laplace transform reads 
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Gr(u!,U2,p)= F[u a =-r,u b =r](ui,u 2 ,p) + F [Ua= _ r ^ Ub=r] (u 1 ,-u 2 ,p) (271) 
K(—T, min(ui, u 2 ))K(max(ui, u 2 ), T) K(—T, —u 2 )K(u\,T) 



K(-T,T) K(-T,T) 
with the function K given by K(u, v,p) = sinh((w — u)^fp) ( [243|) . Here we only need the case u± = = u 2 



(272) 



Gr(0,0,p)=^^) (273) 



This expression may also be found by solving the RSRG equation 

).0:v) = Gl(0.0:v, . . , 

sinh T^fp 



a r G r (0, 0;p) = G r (0, 0;p)B r (0, r,p)B r (T, 0,p)G r (0, 0;p) = G 2 r (0, (hp)—?— (274) 



In the limit of T — > oo, the decomposition (|268| ) becomes 

Z(u L ,u R ;p)= / duE^ (u L ,u;p)E^ (u,u R ;p) + E^(u L , 0;p)G oo (0, 0;p)£;~(0, u R ;p) (275) 



rmin(u L ,u R ) _ i _ 

/ ^ e -(«i-«)VP e -(«K-«)v5 + e- UI -^ — e -WF (276) 



which coincides with (265) as it should. 
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APPENDIX A: EXPLICIT CONSTRUCTION OF THE MEASURE FOR RENORMALIZED LANDSCAPE 

This Appendix is dedicated to the precise definitions and properties of the joint probabilities of extrema introduced 
in Section [n]. 

1. Initial landscape 



The joint probabilities defined in Section II A are precisely defined as follows in terms of the initial landscape 
measure (G4l) 



2n-l 



n AIL J_ 

{x L ,U L ]U 1 ,Xi;...U 2 n-l,X2n-i;X R ,UB)= DUV[U] Y\ S(U Xj -Uj) (Al) 

J 3 = 1 

X\ — 1 X2 — 1 %3 — 1 N—l 

11 8(Ui - U l+1 ) J] 9(U i+1 - Ui) [] 6(Ui - U i+1 ).. Yl 9(U l+1 - U t ) (A2) 



The normalization (26) can be checked by inserting the identity 

N-i 

i = n (^(^+1 - u ^ + 9 ( u * - ( as ) 

i=l 

K xi— 1 X2 — 1 £3 — 1 N—l 

=yi yi n Ui +^ n - u *) n - n - u *) 

n=l l<xi<X2<X3<..<X2n-2<X2„-i<N i=l i=x% i=x 2 i=x 2n -i 



(with K = (N + l)/2, N is odd) into the initial measure 



2. Renormalized landscape 



To define explicitly the measure of renormalized landscapes defined in Section II B , it is convenient to introduce 
the following notation for an ascending bond with x\ < x 2 

X2 — 1 X2 

Qr(xi,x 2 ;U) = J] JJ eir-Ui + UjWUj-UzjBiU^-Ui) (A4) 

i=x± j=i+l 

when U X2 — U Xl > T and otherwise. It expresses that U Xl is the minimum potential over the segment [xi, X2], U X2 is 
the maximum and all descending barriers in the segment with i < j satisfy [/, — Uj < T. Similarly, for a descending 
bond one defines 

X2 — 1 X 2 

Q+( X1>X2 ;U)= J] II OF - U, + Ui)9(U Xl - U 3 )9(U t - U X2 ) (A5) 

i— X\ J— 2+1 

when {/-cj — [/ l2 > r and otherwise. For r = we recover the extrema of the initial landscape 

X2 — 1 

e^ =0 { Xl ,x 2 ;U) = 11 9{U i+l - Ui) (A6) 

x 2 — l 

e+ =0 (x 1; x 2 ; U) = \{ 6{Ui - U l+1 ) (A7) 
As a technicality, to treat the edges we must also introduce 
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2=2 — 1 X-2 

Q+(x u x 2 ;U) = H J] *(T - Uj + WW - U X2 ) (A8) 

i—x\ j=i-\-l 

X-2 — 1 X2 

Q^(x 1 ,x 2 ;U) = J] fl *(r - tf* + 17,- Wi - ( A9 ) 

i— j—i+l 

with no other constraints, since outside the edge sites the potential is at +oo. We may now express the probabilities 
of renormalizcd landscape in terms of the initial landscape as 

~ 2ra-l 

Nr(x L ,UL;u 1 ,Xi;...U2n-i,X2n-i;x R ,UR)= DUV[U] Y[ S{U x .~Uj) (A10) 

J 0=1 

e+(l,x 1 ;U)e^(xi,x 2 ;U)e+(x2,X3}U)..B+(x2n-2,x 2n -i;U)e^(x 2 „-i,N-,U) (All) 



The normalization ( |27| ) can be checked by using the following identity for any F 

K 

l = J2 12 e+{l 7 x 1 :U)Q^(x u x 2 ;U)e+{x 2 ,x 3 ;U)..B+(x 2 n-2,x 2n ^ 1 -U)e^(x 2n ^ 1 ,N-,U) (A12) 

n—l l<Xi<X2-.<X2 n -l<N 

When r is varied, the following important property of renormalization 

X 2 — l X 2 — l 

drQr(xi,x 2 ;U)= £ £ Q-(x 1} yv, U)@+( yi , y 2 ; U)5{U m - U V2 - T)Q^ (y 2 ,x 2 ; U) (A13) 

y 1 =x 1 +l V2=Vl+l 

leads to the renormalization of the full measure (|2 



3. Block product measures 



For the case of block product measure (|3C|), the explicit expressions for the bond measures and for the right edge 
read 

/X3 X 3 — l 

Y[ dU k 6(U X2 -u 2 )S(U X3 -u 3 )0±(x 2 ,x 3 ;U)exp(- ]T Si[U i+1 ,Ui\) (A14) 

k—X2 i—X2 

/X\ Xl — 1 

HdUkSiU^ - Ul )e+(l,x i; C/)exp(- £ 5i[CTi+i.^])*o(Ui) (A15) 
fe=i »=i 

with a symmetric expression for the left edge. 

APPENDIX B: RENORMALIZATION PROCEDURE FOR PERIODIC LANDSCAPE 

In the Section |n|, we have described the renormalization procedure for an interval (xl,xr) with fixed boundary 
conditions. In this Appendix, we briefly describe the changes that are necessary for the case of a periodic landscape 
U\,...Un with Un+i — U\. The number of extrema is necessarily even in the initial landscape as well as in the 
renormalized landscape. The renormalizcd landscape at scale T can thus contain 2n extrema at positions 1 < x± < 
x 2 < .. < x 2n < N and X\ can be either a minimum or a maximum. As V increases there are also samples such that 
all barriers in the systems are smaller than T, corresponding to n = 0. 



The normalization identity analogous to (A12) which allows to treat the periodic case reads 

K 

l = G a {T,U) + Y; Yl {e+(x 1 ,x 2 ;U)Q^{x 2 ,x 3 ;U)..Q+{x 2n - 1 ,x 2n ;U)e^(x2n,xi;U)+ (Bl) 

n=l l<o-i<a;2<..<X2„<A r 

6 r T(xi, x 2 ; U)0^(x 2 , x 3 ; U)..Q^(x 2n ^i,x 2n ; U)0^(x 2n ,xi; U)) (B2) 
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where 

N N 



Go(T, U) = \{ \{ 6{T - Uj + U t )9(T - U t + Uj) (B3) 

i=0 j=0 



is non zero if all barriers in the systems are smaller than T. 

Integrating over the measure J DUV(U), normalized by Z, each term yields: 



K . 

Z = N^ 0) +J2 (N^ n '- ) (u 1 ,x 1 ;...u 2n ,x 2 n) + N^ n ' +) (u 1 ,x 1 ;...U2n,X2n)) (B4) 

1 >/1il...1l9^ 1 ^ _ , ^„ ^ AT 



n=1 •/«!,. .u 2n l<£ci<x 2 <..<a:2„<iV 



where the symbols ± denote whether the bond [x2n>3Ci] is ±- RSRG recursion relations analogous to (|28[ ) may then 
be written, the only new RSRG equation being the growing of flf upon the decimation of the samples with only 
two extrema left. 

Block product measure, involving only Bp can be also written in an obvious way as in (30), and the above 
normalization identity can be written: 

K K 

Z = + Tr(B+B^) n + ^ Tr{B^B+) n (B5) 

n—l n—1 

each term, all integrations excluded, giving the desired TVp 2 ™'^. 

Note that although the starting landscape here can be the same as the one considered in Section |l| by imposing 
U\ = Un in (p5|), the definition of the renormalized landscape here is different from the one in Section [n] for fixed 
boundary conditions. Here one allows consideration of a bond [x2 n ,xi], which can be decimated, whereas in the case 
with infinite left and right boundaries, the decimation of the leftmost and rightmost bonds is not possible. In the 
periodic case, all the bonds can be decimated, resulting eventually in the Go (I 1 ) factor. 



APPENDIX C: DERIVATION OF THE EXPLICIT SOLUTIONS FOR A SCHRODINGER MEASURE 



This Appendix contains the technical details of Section (III) 



1. Solution by a path-integral method 

a. Basic Path-integral 

We first consider the following path integral 

F [UatUh] (u ,x ,u,x) = _ DUex P (- J dx(-( — ) 2 + V[U]) (CI) 

where the sum is over paths which remain in the interval ]u ai Ub[ for xq < x < x\. Seen as a function of (u, x), it 
satisfies the Schrodinger equation 

d x F - d 2 u F - V{u)F (C2) 

with the initial condition at x — x a 

F(u, xq) = S(u — uq) (C3) 

and the absorbing boundary conditions at u — u a and u = Ub 

F(u a ,x)=0 (C4) 
F(u b ,x) = (C5) 



3G 



The Laplace transform with respect to {x — xq) 



F(u,p)= / dxe- p(x - Xo) F{u,x) (C6) 

Jx 



satisfies the system 



d 2 u F - V{u)F =pF- 5(u - uo) (C7) 

F(u a ,p)=0 (C8) 

F(u b ,p)=0 (C9) 

In terms of ((5^), we construct two solutions of the equation (|54| ) which vanish respectively at u a and u b as: 

$-(u,p;u a ) = K(u a ,u,p) (CIO) 

$ + {u,p;u b ) = K(u,u b ,p) (Cll) 

Their wronskian (also independent of u) read: 

W(p;u a ,u b )= $'_{uQ,p;u a )$ + (u Q ,p;u b ) - $_(w ,p; u a )$' + (uo,p; u b ) = K(u a ,u bl p) (C12) 



The solution of the above equations (C9) now read: 

A / | x 1 * / \ / s K{u ai u,p)K(u Q ,u bl p) 

F[u a ,u b ](u,p\u ) = — -<P-{u,p]u a )<P + (u ,p;u b ) = — if u a < u < u (C13) 

W(p;u ai u b ) K{u a ,u b ,p) 

a / i \ 1 ^ , / x K(u a ,u ,p)K(u,u b ,p) 

F\ Ua , Ub \{u,p\u ) = — -$ + (M,p;M b )$_(uo,P;Ua) = 777 \ it U < U < U b (C14) 

W(p;u a ,u b ) K(u a ,u b ,p) 

The probability that the path has remained in the interval ]u a , u b [ at time x 

S[u a ,u„]{x\u ,x Q ) = I duF [UatUb] (u 7 x\u ,x ) (C15) 



satisfies 

- d x S [uaiUb] (x\u ,x ) = rf Ua Ub] (x\u ,x ) + r b [Ua Ub] (x\u ,x ) (C16) 

where and A Ub ^ represent the probability to be absorbed between "time" x and x + dx at the boundary u = u a 
and u = u b respectively. They read: 

r K,«i,] {x\uq,x ) = [d u F [UayUb] (u,x\u ,xo)] \ u=Ua (C17) 

r {u a ,u b ]( x \ u o, x o) = - [0 u -F[„ o ,„„](M,a:K,#o)] |u=„ 6 (C18) 
In Laplace transform with respect to (x — Xq), they have the simple form 

i \ $+(u ,p;u b ) 

v....]fe"°> = Mi^rj < C19 > 



f--i<""»> " 5rS^t <C20) 

We are now in position to evaluate the constrained path-integral (ui, xi;u2,X2) (|50|). 



6. Expression of B T (ui,xi\ U2,X2) for the case U2 — iti + T 

We first consider the function Bp (ux, Xi; U2, X2) for an ascending bond with u% = u\ + T. Setting iti = w Q , the 
distribution of X2 is the distribution of the exit time at the boundary u b = u a + L starting from uq = u a + e without 
being absorbed by the boundary at u a . We thus have to consider the exit probability r? u u ^{x\uq,xq) computed in 
(p20|), for the particular case u$ = u a + e, u b = u a + L 
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l^o (~ r L,n a +r](^l"a + e,x )^ = B [ua Ua+r] (x - x ) (C21) 
where the function B^ u u +r ^(x — Xq) is defined by its Laplace transform 

^K,«a+r](P)= d ^of[u a ,u b ](p;uo)\u =u a ,u b =u +r = + r ,p) (C22) 

c. Expression of Bp (iti, xi; «2, £2) /or t/ie case 1*2 > tti + T 

We now consider the function (ui, x\] U2, X2) for an ascending bond with u 2 > u% + T. To take into account the 
constraint of no returns of more than L, we write the recursion equation 



B T (u ,x ;u + Au,x) = / dx' B r (u , x ;u,x')r b Ua=u _ r Ub=u+Au] (x\u,x') (C23) 

i.e. we obtain in Laplace with respect to {x — x$) the differential equation 

d u B^(u ,;u;p) = Bp (u , ; u;p)d Ub r b [UaUb] (p, u)\ Ua=u -r , Ub =u = (u Q , ;u;p)Z~ (p,u - T,u) (C24) 



with 



7 -f \ d U2 K(u 1 ,u 2 ,p) , 

Z (p,u 1 ,u 2 ) = — ( r— (C25) 

K(Ul,U2,P) 



Thus using the initial condition at u\ = uq + T (C22), we get 

Bp(uo,u,p)— B^ (w ,mo + T,p) exp( / du Z~ (p,u - T,u')) (C26) 

Juo+r 

leading to (»2) given in the text. 



d. Edge bonds 

We now consider the function E^(u,x;ur,xr) for the ascending edge bond near the boundary condition U(x = 
x~^) = +00. Given ur, the only constraint is that u < ur is the deepest minimum with no descending return of more 
than r in between. If u = ur, we simply have x = xr, whereas for u < ur we may again write an equation similar 
to ( C23| ) and we finally obtain in Laplace variable with respect to (xr — x) the result (|5^) given in the text. 



2. Solution from the RSRG equations : Factorized Ansatz and Liouville equation 

The factorized Ansatz (^Tj) for Bp leads to the following RG equation 

r-u — r 

9 r lnBp (u;u";p) = / du'(B^ (V; v! + T;p)) 2 (C27) 

J U 

i.e. from the knowledge of B^(u;u + T;p) for all (T, u) one can reconstruct B^(u;u f, ;p) for all (T,u,i/'). This RG 
equation can be integrated if we introduce an auxiliary function $>(u,v) defined by ( |62| ) since it yields 

pu — r pu 

/ du'd 1 d 2 &(u',u' + T)= dvd^iv-T.v) (C28) 
Ju Ju+r 

u" 

= -dr([ du'd 2 <i>(u' -T,u'))- d r <Z>(u,u + T) (C29) 
Ju+r 
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and thus, integrating (C27) over T one gets 

Br(p;u;u")= exp(-$(u, w + T)) exp(- f du'd 2 $(u' - T,u')) 

Ju+r 

= exp(-$(u" - r, u")) cxp(- [ du'di${u', v! + V)) 

Ju"-r 

which is indeed of the form ( |6l| ) with 

A^{u) = exp(- / du'd^iu', v! + T)) 

A*(u») = cM-*(u" ~ r,u»)) AH J f _ T) 

The consistency of ( |62| ) now requires the famous Liouville equation (|63| ) for the function &(ui,u 2 ). 
As is well known its most general local solution is of the form (Liouville 1850) 

\ 1 ( f^jui) - f R {u 2 ) ^ 
)f R "2 1/2 



(C30) 
(C31) 

(C32) 
(C33) 



(C34) 



where Jl{u) and f R (u) are two arbitrary analytic functions. The alternative form (|64| ) can be obtained by noting 
that the Liouville equation (|63|) implies 



di(e-*dfe*) = 
d^e-^e*) = 



(C35) 
(C36) 



and thus there must exist two functions Vl(u) and V R (u) such that e satisfies the two corresponding Schrodinger 



equations 36 1. The connexion between the two representations is made by noting that fh{u) — ip^ '(it) / 'tpj* (u) ■ 



,(1)/ \ I Nl/2 Il(u) 

Tl W = (- W L) ' 7777 — r 

^ 2 )( U ) = (_^)V2 1 



(C37) 
(C38) 



and similarly for ip R and f R . One also has Vl(u) = j^j — f ( y^y ) 2 ) ano - m g enera l a second solution of the 

Schrodinger equation can be obtained from a first one ^\u) as if>£\u) = V^H/ 11 ^^^)) -2 - Finally one finds 
that one can rewrite (C31) as 



B r (p; u; u"y 



(r L (u)f R (u")) i/2 



exp( 



/fl(«") 



- /fl(tt + r) Jf R (u+T) JlWr) -r)-f R 



(f' L (u)f' R (u"))^ 



fL(u"-V) 



f L (u"-v)-f R {u") ^j fL{u) f L -f R (v(B L ) + ry 



(C39) 
(C40) 



where u(B R ) is the inverse function of f R (u), i.e. f R {u(y)) = y and v (Jl) the inverse function of /l(u). 

Finally, assuming that B^(u; u";p) = A^' {u)A^\u") together with delta function initial conditions, yields for the 
edge blocks 



E+(u L ,u,p) 



E T (u,u R ,p) 



A^{u L ) 
A^{u R ) 



(C41) 
(C42) 
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APPENDIX D: CONFIGURATIONS PRESENTING TWO NEARLY DEGENERATE MINIMA 



1. Distribution of the absolute minimum for the toy model on an arbitrary interval 

The absolute minimum of the toy model on a arbitrary interval [xl, xr] with fixed boundary conditions (ul, ur) is 
given by the T — ► oo limit of the renormalized landscape containing only one valley left ( ps| ) 

/min{u L ,u R ) rXR 
du dxN^ , {x L ,u L - u,x,ur — u,x R ) (Dl) 
-oo J xl 

where 

N^(x L ,u L ~ u,x,u R - u,x R ) = 

e 4 x l-J± XL ( UL -u)ft+( UL _ U) Q; x _ x ^ e+ % XR ( UR -u)-£x 3 R ^ (0, Ui? - U,Xr - x) (D2) 

describes the statistical properties of the position x and the energy u of the absolute minimum of the interval given 
the boundary conditions (ul,xl) and (ur, Xr). To obtain the edge blocks E ( |90| ) in the limit r — * oo, we need the 
function K{u\, u 2 ) in the regime u 2 = Ui + T — ► 00, with u\ fixed. Using the asymptotic expressions of Airy functions 
(|Fq), we thus have 

7T 

K(ui,U2) — —Ai(aui + bp)Bi(av,2 + bp) (D3) 
u 2 -^oo a 

and the edge blocks at T = 00 are simply given by 

^(F,0,p)=^-(0,F,p) = A( ^ ) (D4) 
The function (|D2D describing the distribution of the minimum thus reads 



NV>(x L ,U L - U,X,Ur ~ U,X R ) = e ^2^|-^fl)+f^R(«B-«)-f^i(«i-«) 

/" + °° dAl c'^' 1 -^ 1 A ( a ( ML + U\x) r +co dX 2 jx 2 ( XR -x) MajuR -u) + bi\ 2 ) 
J-oo 2ir Ai(bi\i) 2tt Ai(biX 2 ) { ' 



2. Probability to have two nearly degenerate minima 



To study the measure of configurations with two degenerate minima, we need to consider (|D^) for the special case 
where the absolute minimum u coincides with the left boundary energy Ui 

N^\x L Ax,UR-u L ,XR)^e4(^l) + ^(^) d (x-x L ) [ + °° g» e «V--> \ + 0*) 

For the case ur = Ul + e with e — > 0, this yields for cc^ < x/? 

^(^,0, x, c, x R ) = S(x - x L )e^- x ^ (l + |x R e + 0(e 2 )) 



S(xr — x) + ae 



[ + °° dX 2 lXo(xR _ x) Ai<{UX 2 ) 2 
loo 2tt Ai(WAa) 1 J 



= aed(x - x L )e^ { - Xl x ' r) 



+00 



dX ? ct \2(x R -x L ) Ai, (^2) 

2?r Ai(biX 2 ) 



For the case xr — > 00, after integration over the barrier (it^ -«/,), the formula (Pq) becomes using 



AT^)(a; L ,0,a;,+oo) = lim 

Xr — * + 00 



+oc 



and similarly 



duRN@)(x L ,0,x,UR - u L ,x R ) = e^ XL 5(x - xl)9oo{xl) 
nW(+oo,x,0,xr) = e-^ x ^S(x - x R )g 00 {-x R ) 



(D7) 



(D8) 



(D9) 



(D10) 



Putting together the three results (D8 D9, D10), we get the formula (129) for the configurations with two nearly 
degenerate minima. 
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APPENDIX E: STUDY OF TWO VALLEYS CONFIGURATIONS IN THE TOY MODEL 



This Appendix is dedicated to the study of the two valleys configurations described by (176). Since the barriers F\ 
and F 3 are larger than T satisfying af > 1, we may use the asymptotic expression (F5) for the Bi to get 



N r A) (x u F u x 2 ,F 2 ,x 3 ) 



e ^ 2 ( Fl -F 2 ) f +oo ^ f+oo ^ f +oo ^ f+oo dXi 



r— >oo 7T 
i(\ 1 +\2)x 1 +i(\ 3 -\ 2 )X2-i(\2+>>3) x 3 



27T 



2tt 



2vr 



2tt 



Ai(bi\i)Ai(b\ 2 )Ai(b\ 3 )Ai(bi\i 
At large T, the scaling region for the barrier is thus 

F = r 



(aF 1 + U\' 2 ) 1 /\aF 3 + bl x' 3 )^e-^+ biy ^ /2 -^ aF3+M ^ 



•,3/2 



V 



ay aT 



(El) 
(E2) 

(E3) 



with r\ > of order 1. The barriers are thus concentrated in the vicinity of T. Integration over the two barriers Fi, F 2 
yields at leading order after the rescaling (E3) 



/+00 r+<x> 
dFi J dF 2 Nf ) (xi,F 1 ,x 2 ,F 2 ,x 3 ) 



'r Jr 

■,2{„r\-i/2 



a 2 {aT)- 



-e 3 



{aT) 3 ' 2 _ 



1 



1 



(a 2 X2) 2 
aT 



goc(-xi)goo{xi - x 2 + b(aY) 1/2 )g OQ {x 2 - x 3 + ^aT^^g^xs) 



(E4) 
(E5) 



in terms of (102). 

We now integrate over the maximum x 2 . The constraints of integration x\ < x 2 < x 3 may be forgotten since the 
domains x 2 < x\ and x 2 > x 3 are highly suppressed with the functions g (X] (xi ~ x 2 + b(aT) 1 / 2 ) and #00 (^2 — x 3 + 
fr(ar) 1 / 2 ). The measure for the two minima thus read 



m£\x u x&* ["ax, [ + °°dF 1 [ +e ° dF 2 N?\x u F u x 2 ,F 3 ,x 3 ) * ^(- Xl ) 9oo (x 3 ) e _ |(aF) 3 
J X1 Jr Jr r ^°° a 



/2 



dXi e -iAi(*i+6(«ry/-') f +c° dX2 e -iA 2 (- a3 +6(ar)v-) _ 6(or)1/2|Ai _ A2| 



, 2tt Ai(UXi) 
which yields dl2S|) for aT > 1. 



2?r 



Ai(bi\ 2 ) 



(E6) 
(E7) 



APPENDIX F: USEFUL PROPERTIES OF AIRY FUNCTIONS 



The Airy functions Ai and Bi are two independent solutions of the differential equation 

4f'{z) = z^z) 



of Wronskian 



An immediate consequence is 



AiBi' - BiAi' = - 

7T 



d fBi(x) 



1 



dx \Ai(x) ) TrAi(x) 2 
The asymptotic expressions of the Airy functions at large argument read 



Ai(z) 
Bi(z) 



1 



1 — i _l_3/2 

z 4 e 3 



Re(z)^+oo y>n 



-1/4 ^z 3/2 



(Fl) 
(F2) 

(F3) 

(F4) 
(F5) 
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Another useful property is 



Bi(ix) ,„ , 

lim ; ; ( = ±i F6 

:^±oo Ai(ix) 
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